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8ECTI0N  I 
INTRODUCTION 


Pratt  &  Whitney  Aircraft  Group  Customer  Computer  Deck  (CCD  1144-0.0)  is  a 
three-part  program  consisting  of:  ( 1 )  Sprayring  Model  —  an  analytical  model  that  will  predict 
the  fuel  injection  characteristics  of  turbofan  sprayrings  or  spraybars,  (2)  Flameholder  Combus¬ 
tion  Model  —  an  analytical  model  that  will  predict  the  steady-state  combustion  field  for  a 
turbofan  augmentor,  and  (3)  Rumble  Model  a  dynamic  analytical  model  that  will  predict  the 
rumble  stability  limits  of  turbofan  augmentors.  The  models  may  be  exercised  independently  of 
each  other  or  sequentially  to  supply  data  from  (1)  to  (2)  to  (3).  The  operating  sequence  is 
shown  in  Figure  1. 


FD  223798 


Figure  1.  Combined  Model  Overview 

The  objective  of  the  sprayring  model  is  to  predict  the  fuel  flowrate  distribution  from  a 
specified  fuel  injector  into  the  turbofan  augmentor.  It  is  capable  of  transient  analysis  to 
evaluate  augmentor  throttle  excursions.  The  objective  of  the  flameholder  combustion  model  is 
to  predict  the  augmentor  heat  release  process  for  a  specified  geometry  and  operating  point  in 
a  turbofan  engine.  The  rumble  model  objective  is  to  predict  the  conditions  under  which  the 
specified  turbofan  augmentor  will  experience  low  frequency  instability  (rumble).  Since  the 
dynamic  analysis  is  linear  in  nature,  the  models  predict  rumble  onset  but  not  amplitude. 

The  user  may  execute  any  one  of  the  three  models  independently  of  the  others.  If  the 
sprayring  model  is  not  executed,  the  fuel-air  ratio  distribution  must  be  user  input  to  the 
combustion  model.  If  the  combustion  model  is  not  executed,  then  the  combustion  data  must  be 
user  input  to  the  rumble  model. 


The  User’s  Manual  describes  the  combined  models  and  how  to  use  the  computer  program 
to  predict:  (1)  fuel  system  behavior,  (2)  augmentor  combustion  characteristics,  and  (3) 
turbofan  engine  low  frequency  instability  (rumble)  limits.  To  assist  in  checking  out  the  CCD 
at  the  user  facility,  the  following  items  are  included: 

•  Program  Listings 

•  Test  Case  Input 

•  Test  Case  Output 
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SECTION  II 

TECHNICAL  DISCUSSION 


1.  AUGMENTOR  DESCRIPTION 

Afterburning  is  a  method  by  which  the  maximum  thrust  capability  of  a  basic  engine  may 
be  augmented  by  an  additional  50  percent,  or  more.  Fundamentally,  an  augmentor  (after¬ 
burner)  is  a  ramjet  engine  attached  to  the  turbine  exhaust  case  of  a  turbojet  or  turbofan 
engine.  The  gases  discharged  from  the  turbine  of  the  basic  engine  have  sufficient  velocity  at 
the  higher  thrust  settings  to  satisfy  ramjet  requirements,  regardless  of  whether  the  aircraft  is 
in  a  steep  dive  or  standing  still  at  the  end  of  a  runway. 

The  basic  augmentor  (V-gutter),  Figure  2,  consists  of  only  four  fundamental  parts:  the 
afterburner  duct,  the  fuel  nozzles  or  spraybars,  the  flameholders,  and  a  two-position  or 
variable  area  nozzle.  Because  the  exhaust  nozzle  area  requirements  vary  significantly  depend¬ 
ing  on  whether  or  not  the  augmentor  is  operating,  a  variable  area  exhaust  nozzle  is  in¬ 
corporated. 


Thrust  modulation  in  the  afterburning  mode  is  accomplished  by  varying  the  flow  of  fuel 
to  the  augmentor.  However,  in  order  to  maintain  good  combustion  efficiency  in  the  augmentor 
over  a  wide  range  of  fuel-air  ratios,  the  augmentor  is  separated  into  fuel  supply  “zones”  or 
segments,  for  best  fuel  distribution. 

The  afterburner  duct  must  be  of  such  proportion  that  stable  combustion  can  be 
maintained  during  augmentor  operation.  This  requires  a  burning  section  of  sufficient 
cross-sectional  area  to  ensure  that  the  gas  velocity  through  the  augmentor  does  not  exceed  the 
rate  of  flame  propagation.  Otherwise,  the  flame  would  not  be  able  to  establish  a  firm  foothold 
because  the  onrushing  turbine  exhaust  gases  would  push  the  burning  mixture  right  out  of  the 
exhaust  nozzle.  Fuel  is  introduced  through  a  series  of  perforated  spraybars  located  inside  the 
forward  section  of  the  afterburner  duct.  Not  far  aft  of  these,  flameholders  are  provided  to  help 
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create  local  turbulence  and  to  reduce  the  gas  velocity  in  the  vicinity  of  the  flame.  The 
flameholders  may  take  the  form  of  concentric  rings  or  radial  arms  of  an  angular  “V”  cross 
section,  hence  the  name  V-gutter  augmentor. 

2.  RUMBLE  MODEL  PROGRAM  DESCRIPTION 

a.  General 

The  augmentor  math  model  consists  of  a  set  of  time  dependent  equations  describing  the 
longitudinal  dynamics  of  the  flowing  airstream  and  the  axially  distributed  combustion  process 
in  the  augmentor,  coupled  with  a  solution  technique  for  determining  stability.  These  equations 
are  linearized,  through  the  assumption  of  small  perturbations,  and  transformed  from  the  time 
domain  to  the  Laplace  transform  “S”  domain.  The  solution  technique  is  based  upon  the 
Nyquist  stability  criterion  and  consists  of  determining  whether  the  time  response  of  the 
system  to  a  small  disturbance  would  display  oscillatory  behavior  with  a  growing  amplitude. 
The  result  is  a  determination  of  stability  at  a  given  operating  point,  regions  of  operation  which 
will  cause  rumble,  and  changes  to  the  augmentor  to  make  it  rumble-free. 

b.  Modeling  Approach 

Since  rumble  has  long  been  associated  with  the  relatively  low-frequency  longitudinal,  or 
axial,  modes  of  vibration  of  the  air  column  in  the  augmentor,  the  model  was  formulated  to  take 
only  the  longitudinal  dimension  into  account.  Accordingly,  each  station  in  the  model  was 
considered  to  reprsent  a  plane  over  which  the  value  of  any  parameter  (such  as  velocity, 
pressure  or  density)  could  be  considered  as  uniform  at  any  instant  in  time. 

Motion  pictures  of  rumble  have  shown  a  change  in  color  of  the  burning  gas  during  a  '-ycle 
of  oscillation,  indicating  that  alternate  hotter  and  cooler  combustion  products  were  being 
produced.  These  hot  and  cold  combustion  products  could  be  seen  drifting  from  the 
flameholder  to  the  exhaust  nozzle  in  a  time  span  which  matched,  or  was  a  multiple  of,  the 
period  of  oscillation  of  the  rumble.  Since  flowrate  out  through  the  nozzle  is  dependent  upon 
the  temperature  of  the  entering  gas,  it  was  important  that  the  model  treat  the  traveling 
combustion  products,  which  were  mathematically  identified  as  traveling  entropy  waves. 

The  equations  developed  for  describing  rumble  can  be  classified  into  two  types.  First, 
there  are  the  momentum,  continuity  and  energy  equations,  together  with  the  boundary 
conditions,  which  describe  how  each  parameter  at  any  station  in  the  augmentor  responds  to  a 
disturbance  in  combustion  heat  release.  These  are  referred  to  as  the  acoustic  equations. 
Secondly,  there  are  the  combustion  equations  which  describe  how  combustion  heat  release 
responds  to  variations  in  the  system  parameters  such  as  velocity,  pressure  and  density. 
Together,  the  acoustic  equations  and  the  combustion  equations  describe  the  rumble  mecha¬ 
nism,  by  which  a  disturbance  in  combustion  causes  a  disturbance  in  velocity,  pressure  and 
density  throughout  the  augmentor  which  in  turn  causes  a  disturbance  in  combustion.  A 
description  of  the  equations,  boundary  conditions  and  assumptions  is  presented  in 
Appendix  A. 

Since  the  purpose  of  the  program  was  to  develop  an  understanding  of  the  rumble 
mechanism  and  demonstrate  that  the  onset  of  rumble  could  be  predicted,  thereby  defining  the 
boundary  between  stable  and  unstable  operating  regions,  it  was  necessary  only  to  model  the 
augmentor  for  the  first  few  increments  of  time  before  the  oscillation  had  built  up  into  an 
appreciable  amplitude.  This  allowed  use  of  a  small  perturbation  technique  which  led  to  linear 
equations  and  mathematical  simplification.  Linear  equations  can  describe  the  system  for  small 
oscillation  amplitudes  and  can  predict  whether  the  system  initially  at  rest  would  begin  to 
oscillate.  Because  the  nonlinearities  associated  with  large  amplitude  oscillations  (which  even¬ 
tually  stop  the  amplitude  from  growing)  were  ignored,  the  linear  equations  do  not  allow  a 
prediction  of  the  final  limit -cycle  amplitude. 
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c.  Model  Description 


The  rumble  model  was  designed  for  simple  input-output  and  requires  no  intermediate 
engineering  interpretation,  Figure  3.  The  input  requires  engine  geometry  and  pressures, 
temperatures  and  Mach  numbers,  obtained  from  engine  steady-state  cycle  tables.  The  user 
may  select  to  input  augmentor  fuel-air  ratio  and  empirical  combustion  data  or  he  may  exercise 
the  flameholder  combustion  model  which  calculates  and  supplies  the  required  augmentor 
combustion  data  to  the  rumble  model.  An  input  option  allows  the  user  to  specify  the  specific 
augmentor  types  (V-gutter,  Vorbix  or  Swirl)  to  be  used.  No  calculation  or  dynamic  information 
is  required.  The  user  may  select  either  tabular  and  plotted  output  or  only  plotted  output,  as 
shown  in  Figure  4.  From  the  plot  the  user  identifies  the  frequencies  at  which  the  phase  is  zero. 
He  then  checks  the  gain  at  each  of  the  identified  frequencies.  If  the  gain  is  one  or  greater,  the 
program  has  predicted  that  rumble  will  occur.  If  the  gain  is  less  than  one,  the  program  has 
predicted  that  the  operating  point  is  stable.  For  example,  Figure  4  indicates  rumble  at  60  Hz 
and  at  140  Hz.  The  user  can  then  change  geometry  or  operating  point  inputs  and  repeat  the 
process  to  determine  the  effects  of  the  change.  This  form  of  output  was  chosen  because  it 
facilitated  development  of  the  model,  yielded  a  compact,  easy  to  interpret  answers,  and  made 
better  use  of  computer  time  than  a  time-domain  solution. 

To  model  rumble  required  a  transient  description  of  the  longitudinal  dynamics  of  a 
turbofan  engine.  To  computerize  the  formulation,  the  mathematical  description  was  simplified 
by  restricting  the  range  of  validity  of  the  equations  to  small  perturbations  about  a  mean 
steady-state  operating  point.  This  allowed  linearization  of  the  equations  to  a  form  which 
correctly  describes  small  scale  transients,  but  in  which  the  nonlinear  terms  which  are 
important  in  large  scale  transients  could  be  omitted.  The  resultant  linearized  equations 
describe  the  initial  period  of  the  time  when  rumble  oscillations  begin  to  grow  and  are  valid  to 
the  point  where  the  rumble  amplitude  reaches  values  at  which  the  nonlinear  terms  become 
important.  This  is  sufficient  to  determine  whether  an  engine,  if  placed  at  a  given  operating 
point,  will  spontaneously  bloom  into  rumble. 

The  model  could  have  been  made  to  yield  solutions  in  the  time-domain  by  programming 
the  equations  on  an  analog  computer.  The  output  would  have  been  a  time  trace  of  any  selected 
parameter  (e.g.,  augmentor  pressure).  At  a  stable  operating  point  the  trace  would  be  a  straight 
line,  whereas  at  an  unstable  operating  point  the  trace  would  show  a  sinusoidal  oscillation  with 
an  increasing  amplitude.  The  amplitude  would  grow  without  bound  for  as  long  as  the  solution 
continued,  because  of  the  omission  of  the  nonlinear  terms. 

The  same  information  can  be  more  easily  obtained  by  a  nontime-domain  solution 
technique.  Such  a  technique  was  chosen  for  the  rumble  model.  Commonly  called  the  Nyquist 
criterion,  it  is  based  upon  the  fact  that  the  allowable  forms  of  the  time  domain  solution  are 
known.  This  technique  allows  use  of  a  matrix  program  which  can  quickly  solve  large  numbers 
of  simultaneous  equations. 

The  Nyquist  criterion  is  a  procedure  which  makes  use  of  the  Laplace  transform  and 
conformal  mapping  to  determine  whether  the  transient  solution  would  show  unstable  behavior. 
To  apply  the  criterion,  the  time-domain  equations  are  transformed  into  the  Laplace  “S” 
domain.  The  result  is  a  square  homogenous  matrix.  The  determinant  of  the  matrix  coefficients 
is  a  function  of  “S”,  called  the  characteristic  function,  and  contains  all  of  the  information 
needed  to  determine  whether  the  system  being  described  is  stable  or  unstable.  If  all  zeros  of 
the  characteristic  function  have  negative  real  parts,  the  system  is  stable;  if  any  zeros  have 
positive  real  parts,  the  system  in  unstable.  Conformal  mapping  is  used  to  examine  the 
characteristic  function  for  the  presence  of  zeros  with  positive  real  parts. 
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Figure  3.  Rumble  Model  Flow  Diagram 
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Figure  4.  Rumble  Model  Plotted  Computer  Output 

To  accomplish  the  conformal  mapping  the  equations  which  describe  the  augmentor  were 
written  to  describe  a  “feedback  loop.”  The  feedback  loop  was  formed  for  the  rumble  model  by 
considering  that  the  combustion  rate,  called  q’jn  was  an  “input”  to  the  acoustic  equations.  This 
yielded  as  “output”  the  pressure,  velocity  and  density  at  each  station  throughout  the  engine. 
The  “output”  was  then  considered  to  feedback  through  the  combustion  equations  to  form  a 
“feedback”  combustion  heat  release  rate,  called  q’OUf  The  resultant  “loop”  is  shown  in  Figure 
5.  Actually,  only  one  heat  release  rate  is  present.  The  use  of  the  two  names  q’jn  and  q’out  allows 
the  formation  of  the  ratio  q’„ut/q’in*  called  the  "Open  Loop  Transfer  Function”  (OLTF). 
COnformal  mapping  to  examine  the  zeros  of  the  characteristic  function  is  carried  out  by  using 
the  OLTF. 

Referring  to  Figure  5,  the  heuristic  argument  can  be  made  that  if  a  loop  is  subjected  to 
an  externally  supplied  sinosidal  input  (q’jn)  and  it  returns  a  feedback  (q’out)  which  is  in  phase 
with  the  input  (<£=  0)  and  of  equal  amplitude  (gain  =1),  then  the  externally  supplied  input 
could  be  removed  and  the  loop  would  continue  to  oscillate.  A  gain  greater  than  one  then 
implies  that  the  loop  would  be  driven  to  even  higher  amplitude,  while  a  gain  less  than  one 
implied  that  the  oscillations  would  die  out  once  the  input  was  removed.  The  model  determines 
whether  the  time  solution,  if  calculated  would  display  oscillatory  behavior  with  a  growing 
amplitude.  It  does  this  through  a  solution  technique  which  is  simpler  and  faster  to  apply  than 
a  solution  in  the  time-domain. 
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Figure  5.  “Feedback  Loop"  Visualization  of  Rumble  Model 

3.  FLAMEHOLDER  COMBUSTION  MODEL  PROGRAM  DESCRIPTION 

a.  General 

The  combustion  model  performs  a  multistreamtube  analysis  of  the  flame  stabilization 
and  propagation  phenomena  in  a  turbofan  augmentor.  The  augmentor  is  divided  into  a 
multitude  of  equivalent  two-dimensional  streamtubes  with  a  single  flameholder  element  in 
each.  The  program  evaluates  each  streamtube  and  then  mass  averages  the  results. 

For  each  streamtube  the  program  proceeds  from  the  augmentor  inlet  towards  the  exhaust 
nozzle  and  evaluates  each  step  in  the  stabilization  and  propagation  of  the  augmentor  process. 
The  ultimate  result  is  the  level  of  combustion  efficiency  in  that  streamtube.  The  program  then 
performs  a  small  perturbation  in  velocity,  pressure,  inlet  temperature  and  fuel-air  ratio  to 
evaluate  the  efficiency  slopes. 

The  final  outputs  are  the  fan  duct  efficiency,  the  core  stream  efficiency  and  the  efficiency 
slopes  with  respect  to  the  four  perturbed  variables. 

b.  Modeling  Approach 

The  approach  taken  for  each  streamtube  is  a  step-by-step  solution  to  the  physical 
phenomena  which  determine  the  flame  stability  limits  of  the  spraybar-flameholder  configura¬ 
tion  and  the  subsequent  turbulent  flame  propagation  rate.  These  phenomena  include  liquid 
fuel  injection,  droplet  formation,  vaporization,  fuel  impingement  onto  the  flameholder,  wake 
reaction  kinetics  and  turbulent  flame  penetration. 
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The  approach  used  is  different  for  the  fan  duct  streamtubes  and  the  core  streamtubes. 
The  necessity  for  different  approaches  lies  in  the  degree  of  liquid  fuel  vaporization  between 
the  spraybar  and  the  flameholder.  In  the  core  streamtubes,  the  fuel  is  virtually  totally 
vaporized  in  the  first  few  inches  by  the  hot  turbine  exhaust  flow.  In  the  fan  duct  stream,  the 
much  cooler  airflow  results  in  only  a  slight  degree  of  vaporization  in  the  four  to  six  inches 
typical  spraybar  to  flameholder  distance. 

The  core  stream  analysis  is  thus  done  assuming  that  the  fuel  at  the  flameholder  is  in  the 
vapor  phase  and  the  flameholder  wake  fuel-air  ratio  is  the  same  as  the  total  fuel-air  ratio.  This 
value  is  used  in  the  kinetics  analysis  of  the  wake  reaction  to  evaluate  the  stability  limits. 

In  the  fan  duct  streamtubes,  however,  the  low  level  of  droplet  vaporization  yields  a  vapor 
phase  fuel-air  ratio  at  the  flameholder  which  is  well  below  the  lean  limit  for  hydrocarbon  fuels. 
Since  the  liquid  fu^.  droplets  are  not  capable  of  entering  the  flameholder  recirculation  wake 
due  to  their  excessive  momentum,  there  must  be  some  other  mechanism  to  provide  the 
necessary  wake  vapor  fuel  for  stable  combustion. 

This  mechanism  in  the  fan  duct  streamtubes  is  the  collection  of  the  liquid  fuel  droplets 
onto  the  surface  of  the  flameholders  and  the  vaporization  of  the  resultant  liquid  film.  This 
evolved  vapor  recirculates  into  the  flameholder  wake  with  a  portion  of  the  droplet  evolved 
vapor  fuel  to  generate  the  wake  vapor  fuel  concentration. 

The  streamtube  analyses  compute  the  degree  of  wake  reaction  at  the  level  of  vapor 
fuel-air  ratio  appropriate  to  the  streamtube  type  and  approach  conditions.  For  the  fan  duct 
cases,  this  requires  a  convergent  solution  between  the  wake  kinetics  and  the  surface  vaporiza¬ 
tion. 


Once  the  flameholder  wake  reaction  is  evaluated,  the  analysis  computes  the  rate  of  flame 
penetration  into  the  free-stream  as  a  turbulent  flame  sheet.  This  rate  is  adjusted  by  the  wake 
reaction  level  to  account  for  the  ignition  response  in  the  recirculation  zone  shear  layers.  The 
flame  penetration  rate  is  integrated  over  the  available  augmentor  length  to  provide  the  level  of 
streamtube  efficiency. 

The  program  thus  performs  a  quantitative  evaluation  of  the  phenomological  processes 
which  occur  in  the  turbofan  augmentor.  The  individual  calculations  are  a  combination  of 
analytical  evaluations  and  empirical  results  as  required  to  ensure  quantitative  accuracy. 

c.  Model  Description 

The  combustion  model  was  designed  as  a  complete  unit.  The  program  does  not  require 
on-line  engineering  interaction.  The  combustion  model  may  be  run  as  a  separate  entity  or  as 
a  generator  for  subsequent  stability  analysis  with  the  rumble  model.  When  exercised  alone,  the 
combustion  model  is  an  augmentor  analysis  program  and  the  output  is  a  comprehensive 
description  of  the  injection,  stabilization  and  flame  propagation  processes.  In  this  mode,  the 
program  is  useful  as  a  design  tool  for  conventional  turbofan  augmentors.  The  effects  of  fuel 
system  distribution  and  V-gutter  flameholder  tailoring  may  be  determined. 

When  exercised  in  conjunction  with  the  stability  analysis,  a  less  extensive  output  is  given 
and  the  prime  purpose  of  the  program  is  to  generate  the  response  of  augmentor  efficiency  to 
variations  in  fuel-air  ratio  and  inlet  velocity,  pressure  and  temperature. 

The  augmentor  breakdown  and  specific  description  of  one-stream-tube  is  shown  in  Figure 
6.  For  a  single  fan  duct  streamtube,  the  computer  logic  is  shown  in  Figure  7.  The  identified 
subroutines  each  evaluate  a  specific  portion  of  the  overall  combustion  process. 
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Figure  6.  Location  of  a  Core  Streamtube  in  a  Turbofan  Engine  Augmentor 
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Figure  7.  Single  Streamtube  Logic  Map 
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The  input  requirements  for  a  fan  duct  streamtube  are  those  to  fully  describe  the 
approach  flow  field,  geometry  of  the  streamtube  and  ftameholder,  and  the  total  fuel-air  ratio. 
The  execution  of  one  streamtube  proceeds  as  follows: 

(1)  Inject 

This  subroutine  evaluates  the  droplet  sizes  formed  by  a  variable  area  spraybar  as  a 
function  of  the  injection  pressure  drop.  Five  droplet  sizes  are  calculated  which  represent  the 
cumulative  volume  versus  pressure  drop  curve  for  this  spraybar  type. 

This  subroutine  evaluates  the  amount  of  the  liquid  fuel  which  is  flash -vaporized  by  the 
injection  process.  This  evaluation  is  performed  as  an  adiabatic  expansion  process  from  the 
high-pressure  spraybar  fuel  condition  to  the  low-pressure  augmentor  conditions.  The  ap¬ 
propriate  fuel  enthalpy  chart  is  used,  keyed  by  the  fuel  type  input  variable. 

The  liquid  flowrate  which  remains  is  partitioned  equally  into  the  five  size  groups.  The 
total  flowrate  is  originally  calculated  from  the  total  fuel-air  ratio  input  and  the  airflow  which 
is  calculated  from  the  streamtube  geometry  and  flow  conditions. 

(2)  Accel 

This  program  subroutine  evaluates  the  rate  of  droplet  vaporization  and  acceleration 
which  occurs  between  the  spraybar  and  the  downstream  V-gutter  flameholder. 

The  equations  for  acceleration  assume  a  spherical  liquid  droplet  which  is  accelerated  by 
drag  forces  only.  The  drag  coefficient  is  evaluated  as  a  function  of  Reynold’s  number  based  on 
the  relative  air-liquid  velocity. 

Concurrently,  the  rate  of  liquid  vaporization  is  evaluated  as  forced  convection  mass 
transfer  utilizing  a  mass  transfer  Nusselt  number  correlation  which  is  also  based  on  the 
relative  velocity  Reynold’s  number.  The  requirement  to  simultaneously  solve  the  vaporization 
and  acceleration  equations  was  met  by  a  finite  difference  solution.  A  small  time  increment  is 
selected  and  the  acceleration  solution  performed  to  generate  a  velocity  increase  for  the  liquid 
droplet.  Using  the  average  velocity  over  this  time  increment,  a  vaporization  rate  is  calculated 
and  a  vaporized  fraction  evaluated.  This  sets  a  new  droplet  size  for  the  next  time  interval.  The 
average  velocity  over  this  time  is  also  used  to  calculate  a  distance  travelled. 

This  procedure  is  repeated  until  either  the  liquid  droplet  reaches  the  flameholder  or  is 
fully  vaporized.  This  analysis  is  repeated  for  each  size  group  of  the  five  initially  set. 

(3)  Collect 

At  the  flameholder  plane,  the  program  evaluates  the  rate  of  liquid  deposition  onto  the 
surface  of  the  V-gutter.  This  deposition  occurs  as  the  liquid  droplets  are  unable  to  follow  the 
divergent  air  flow  streamlines  around  the  leading  edge  of  the  flameholder. 

The  evaluation  of  the  rate  of  deposition  is  performed  as  a  correlative  solution  to  the  point 
where  liquid  droplets  just  hit  the  flameholder  surface.  The  variables  include  flameholder 
geometry,  droplet  diameter  and  flow  conditions.  The  correlation  equations  are  based  on 
calculations  which  were  done  externally  to  this  program,  where  limit  trajectories  were  estab¬ 
lished  based  on  potential  flow  solutions  to  the  flow  field  approaching  the  flameholder. 
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The  program  utilizes  the  droplet  diameter  which  exists  after  the  vaporization  evaluation 
to  calculate  the  percentage  of  the  liquid  flowrate  in  each  size  group  which  is  deposited  on  the 
V-gutter  surface.  This  is  done  for  each  of  the  five  size  groups.  The  collection  mass  flowrate  is 
evaluated  from  each  size  group  collection  percentage  and  the  liquid  flowrate  in  each  group  at 
the  flameholder. 

(4)  Recirc 

The  gaseous  recirculation  rate  into  the  flame  holder  wake  is  evaluated  from  a  variety  of 
literature  sources  which  present  recirculation  zone  volume  and  flowrate  as  a  function  of 
flameholder  geometry  and  flow  conditions.  The  program  evaluates  a  “recirculation  efficiency” 
which  is  the  ratio  of  recirculated  mass  flow  to  the  flowrate  through  the  area  blocked  by  the 
flameholder.  This  typically  runs  15  to  25%. 

The  correlations  cover  a  range  of  the  variables  which  control  the  recirculation  such  as 
flameholder  apex  angle,  blockage,  approach  Mach  number,  and  temperature.  The  result  of  the 
subroutine  is  the  recirculation  zone.  These  are  used  in  the  analysis  of  the  wake  reaction 
efficiency. 

(5)  Wake 

The  wake  reaction  is  treated  as  if  it  occurred  in  a  well-stirred  reactor  with  volume  and 
entry  flowrate  as  evaluated  in  RECIRC.  The  kinetics  are  assumed  to  proceed  as  a  single-step, 
second  order  conversion  process.  The  kinetics  utilize  rate  coefficients  which  simulate  aircraft 
fuel  behavior.  The  required  inputs  are  wake  volume,  wake  fuel-air  ratio,  recirculation  rate  and 
inlet  conditions  of  pressure,  temperature,  etc.  The  output  of  the  analysis  is  the  wake  reaction 
efficiency  and  mean  wake  temperature. 

(6)  Beta  3 

This  subroutine  evaluates  the  degree  of  vaporization  of  the  liquid  film  which  exists  on  the 
flameholder  surface.  The  vaporization  process  is  one  forced  convection  from  the  surface  into 
the  trailing  wake  shear  layer  and  heat  transfer  from  the  flameholder  wake  through  the 
flameholder  metal  into  the  liquid  film.  The  program  utilizes  a  small  element  approach  using  10 
elements  on  each  side  of  the  flameholder.  The  mass  flux  and  heat  flux  are  evaluated  for 
one-at-a-time  starting  at  the  flameholder  leading  edge.  Any  liquid  left  unvaporized  is  assumed 
to  leave  the  trailing  edge  of  the  flameholder  and  traverse  through  the  wake  shear  layers 
downstream. 

The  solution  of  WAKE  and  BETA  3  must  be  done  simultaneously  since  BETA  3 
requires  wake  temperature  to  find  fuel  vaporization  and  the  vaporization  influences  WAKE 
through  fuel-air  ratio.  The  solution  approach  is  described  in  Appendix  B  with  a  typical  result 
shown  here  in  Figure  8. 

(7)  Flame 

The  turbulent  flame  propagation  downstream  of  the  flameholder  uses  a  small  step 
difference  solution  with  axial  profiles  of  turbulence,  flow,  etc.  The  procedure  is  described  in 
Appendix  B. 
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Figure  8.  Duct  Stream  Flameholder  Wake  Solution 

d.  Input  Requirements  and  Comments 

The  model  requires  as  input  the  physical  variables  which  describe  the  fan  duct  and  core 
stream  geometry  and  operating  conditions.  Since  the  model  functions  by  repetitive  analysis  of 
single  streamtubes,  the  input  is  required  for  each  different  type  of  streamtube.  A  different 
type  is  one  with  any  input  variable  different. 

The  input  requires  the  following  values  along  with  the  input  described  in  Section  11.7. 

BPR  Actual  value.  Default  to  1.0  if  run  as  a  duct  burner  with  no  core 

engine  and  WCOOL=l/2  X  (mcool/mduct) 

M6C  Inlet  Mach  numbers 

MfiH 

NTC  No.  of  types  of  fan  duct  streamtubes 

NTH  No.  of  types  of  core  streamtubes 

PS6  Inlet  static  pressure,  psia 

Array  input  is  required  to  describe  each  streamtube  fully.  These  array  values  are 
aerodynamic  and  geometric.  The  array  is  the  number  of  streamtube  types  in  the  fan  (NTC)  or 
core  (NTH)  and  the  number  of  stream  flow  tubes  of  each  type  identified  in  the  fan  (NSC)  or 
core  (NSH)  sections.  If  three  different  types  of  fan  streamtubes  are  used  (NTC  =  3),  with  a 
total  number  of  28  fan  streamtubes  (18  of  flameholder  width  (FHWC)  =  1.0  in.,  4  of 
flameholder  width  =  0.75  in.  and  6  of  flameholder  width  =  1.25  in.,  NSC  =  18,  4,  6)  and  if  the 
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first  two  types  operate  at  the  same  fuel-air  ratio  (FAC),  but  different  from  the  third,  then  the 
input  to  the  model  to  describe  this  case  would  be: 

dlnput . ,  NTC  =  3,  NSC  =  18,  4,  6,  FHWC  =  1.0,  0.75,  1.25,  FAC  =  0.05, 

0.05,  0.045,.... 

SEnd 

The  program  as  currently  written  assumes  a  unit  depth  streamtube,  i.e.,  1  inch  depth. 
The  mass  flowrates  will  be  based  on  this  value.  If  true  flow  values  are  required,  the  number  of 
each  type  (NSC  or  NSH)  should  be  the  number  of  1  inch  deep  streamtubes  of  that  type.  This 
is  shown  in  Figure  9. 
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^—Typical  f"an  Duct  Streamtube 


Figure  9.  Location  of  Typical  Fan  Duct  Streamtube 

The  geometric  inputs  required  for  a  single  streamtube  are  shown  in  Figure  10.  The  value 
of  blockage  is  referenced  to  Figure  9.  The  input  should  reflect  the  ratio  of  flameholder  width 
to  the  streamtube  limits.  This  value  of  blockage  sets  the  required  flame  penetration  for  100% 
efficiency  and  must  be  input  correctly. 

The  value  of  EPSC  is  the  approach  turbulence  and  will  affect  the  flame  speed.  Unless 
specific  data  are  available,  use  a  value  of  0.04  for  a  turbofan  engine. 

The  input  value  for  PFSR  controls  the  mean  droplet  size  from  the  sprayring,  which  has 
data  from  a  variable  area  orifice  built  in.  If  other  values  are  desired,  use  the  equation: 

da,  =  795  —  (PFSR  —  PS6)  -4 
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to  determine  the  input  value  of  PFSR  required  to  yield  a  desired  mean  droplet  diameter,  in 
microns.  This  is  the  only  place  where  PFSR  is  used,  so  no  disruption  occurs  if  nontrue  values 
are  wiput. 


tO  >M4M 

Figure  10.  Single  Streamtube  Geometry  and  Flow  Inputs  —  Fan 

For  the  aerodynamic  inputs,  also  reference  Figure  10,  the  required  input  is  shown.  As 
previously  mentioned.  PS6  is  assumed  to  be  uniform  across  the  streamlubes. 

One  input  set  requires  external  evaluation.  This  is  the  values  assigned  to  WEXT  and 
TEXT  in  the  fan  duct  streamtubes.  The  purpose  of  this  input  is  to  account  for  the  influence 
of  hot  gas  migration  down  the  wake  region  of  the  fan  duct  flameholders  from  either  the  core 
or  from  a  pilot.  WEXT  is  defined  as  the  ratio  of  this  “external”  flowrate  to  the  recirculated 
flowrate.  To  allow  for  flexibility  in  design  selection,  this  input  format  was  selected.  The  user 
must  evaluate  whatever  flowrate  is  expected  and  calculate  WEXT.  For  use  in  estimating  the 
recirculated  flowrate,  assume  K,  =  0.25  use: 

mr  =  K|  -pa  VaN 

for  recirculation  rate  per  inch  of  flameholder  length.  Typical  values  of  WEXT  are  0.02  to  0.04. 
TEXT  is  the  temperature  of  this  "external”  flowrate.  These  are  shown  in  Figure  11. 

The  liner  cooling  airflow  input,  WCOOL,  is  the  ratio  of  cooling  liner  air  to  total  engine 
air.  As  such,  the  engine  bypass  ratio  is  required  to  evaluate  the  net  available  fan  duct  airflow. 
If  no  cooling  air  is  taken  from  the  fan  duct  or  if  input  fuel-air  ratios  are  based  on  the  true  net 
air  available  for  combustion,  input  WCOOL  =  0.0.  If  a  duct  burner  is  being  analyzed  and  it 
does  have  a  cooling  liner,  set  a  dummy  valve  of  1.0  for  BPR  and  set  WCOOL  by: 

WCOOL  -  -cooling/,  m„ucl 


Core  Flameholder 
Gutter 


e.  Output 

The  program  has  two  output  formats,  long  and  short.  The  long  format  presents  detailed 
values  for  the  processes  which  control  the  wake  vapor-phase  fuel-air  ratio  and  flame  penetra¬ 
tion.  The  short  format  essentially  presents  the  overall  results.  For  both,  the  results  are 
presented  as  a  streamline-by-streamline  analysis  with  fan  and  core  summaries. 

(1)  Fan  Streamtube 

The  long  format  presents  the  input  data  for  each  streamtube  and  two  calculated  values. 
These  values  are  the  effective  streamtube  fuel-air  ratio  and  the  effective  recirculation 
temperature.  The  equations  used  for  these  are  shown  in  Appendix  B. 

The  output  lists  the  calculated  values  for  the  injection  process,  mean  droplet  size  and 
flash  vaporization,  and  the  influence  coefficients,  di,  $■><  d:i  and  K|,  which  control  the  wake 
fuel-air  ratio.  The  importance  of  these  values  is  explained  earlier  in  this  section  and  in 
Appendix  B. 

A  word  of  caution  is  in  order.  If  the  output  is  preceded  by  the  warning  that  the  wake 
temperature  iteration  has  failed,  the  situation  is  such  that  the  wake  has  exceeded  the  rich  limit 
at  the  input  conditions.  Although  output  is  presented,  it  is  not  valid  and  merely  represents  the 
limits  of  the  internal  convergence  search  routine.  For  example,  wake  temperature  will  always 
be  5(XX)0F  for  a  failed  case.  If  a  single  streamtube  is  being  run,  several  other  error  messages 
will  result  as  the  program  attempts  to  interpret  zero  efficiency.  If  multiple  streamtubes  are 
being  run,  the  program  will  ignore  the  failed  streamtube  in  all  calculations. 

The  initial  flame  speed  is  the  laminar  flame  speed  at  the  appropriate  inlet  conditions. 
The  turbulence  level  is  the  value  induced  by  the  flameholder. 


In  the  stream  efficiency  section  for  each  streamtube,  the  following  comments  are 
applicable: 

•  The  ideal  temperature  use  is  based  on  the  effective  fuel-air  ratio. 

•  The  efficiency  is  the  ratio  of  flame  penetration  to  streamtube  width  at 
the  exhaust  nozzle. 

•  The  actual  temperature  rise  is  based  on  the  above  conditions.  The  exit 
temperature  is  based  on  streamtube  inlet  plus  this  actual  temperature 
rise. 

•  The  flowrates  are  for  a  1  inch  deep  streamtube.  The  fuel  flowrate  uses  the 
effective  fuel -air  ratio. 

The  fan  streamtube  summary  presents  the  major  items  from  each  streamtube  and  then 
the  exit  average  results.  The  cooling  air  flowrate  ratio  is  repeated  here.  Two  more  values  of 
combustion  efficiency  are  presented  and  two  values  of  average  exit  temperature. 

The  average  streamline  exit  temperature  is  the  mass  weighted  average  of  the  individual 
exit  temperatures.  The  chemical  combustion  efficiency  is  based  on  this  value  for  exit  and  an 
ideal  temperature  use  based  on  the  average  effective  fuel-air  ratio  and  average  inlet 
temperature. 

The  average  duct  exit  temperature  includes  the  mass  weighted  effect  of  the  liner  cooling 
air  being  added  to  the  streamtubes  at  the  exhaust  nozzle  inlet.  The  average  thermal  combus¬ 
tion  efficiency  is  based  on  this  exit  temperature,  the  average  inlet  temperature  and  an  ideal 
temperature  rise  is  based  on  the  average  input  fuel-air  ratio. 

Since  engine  analysis  procedures  generally  base  the  fan  duct  fuel-air  on  the  total  duct 
airflow  and  use  the  thermal  nozzle  inlet  averages,  the  value  of  thermal  combustion  efficiency 
is  the  one  which  is  used  for  rumble  prediction. 

The  total  flowrate  presented  here  includes  the  number  of  each  type  of  streamtube  as  do 
all  of  the  above-mentioned  mass  averaged  values. 

Also  note  that  at  no  time  are  efficiencies  ever  mass  averaged  directly.  All  average 
efficiencies  are  based  on  comparison  of  the  average  results  of  individual  streamtubes  to  the 
result  of  the  average  inlet.  That  is,  the  burn-then-mix  process  is  compared  to  the  ideal 
mix-then-burn  process.  Since  curves  of  ideal  temperature  rise  exhibit  peak  vs.  fuel-air  ratio, 
the  average  efficiency  of  two  streamtubes,  one  lean  and  one  rich,  may  very  well  be  less  than 
either  streamtube  separately. 

(2)  Core  Streamtubes 

Due  to  the  absence  of  droplet  effects,  the  output  is  greatly  sim;  lified.  The  wake  reaction 
results  are  presented  as  well  as  initial  flame  properties.  Without  liner  cooling  air,  there  is  no 
fuel-air  ratio  shift  and  thus,  only  one  efficiency  definition.  The  process  of  evaluation  of  the 
ideal  temperature  rise  is  given  in  Appendix  B.  All  of  the  comments  in  the  fan  stream  apply 
here,  except  that  thermal  efficiency  is  not  defined  here  due  to  the  lack  of  liner  cooling  air. 

If  the  message  “Aerodynamic  Loading  exceeds  Kinetic  Capacity"  occurs,  the  blowout 
limits  were  exceeded. 


4.  SPRAYRING  MODEL  PROGRAM  DESCRIPTION 

a.  General 

The  sprayring  model  is  a  transient  analysis  which  solves  the  hydrodynamic  and 
thermodynamic  equations  for  fuel  tlow  into  and  out  of  an  augmentor  fuel  injection  system. 
The  type  of  fuel  system  may  be  either  spraybars  or  sprayrings.  Various  options  are  available 
to  specify  the  manner  in  which  the  fuel  flowrate  into  the  injection  system  is  controlled.  The 
model  is  capable  of  evaluating  up  to  sixteen  separate  zones  of  fuel  injection.  The  output 
consists  of  the  fuel  flowrate  distribution  at  a  specified  transient  time.  The  output  is  used  to 
define  the  fuel-air  ratio  distribution  for  later  use  by  the  flameholder  combustion  model.  The 
ouput  may  also  be  used  directly  to  evaluate  fuel  injector  designs  by  comparison  of  transient 
pressures  and  fuel  distributions. 

b.  Modeling  Approach 

The  augmentor  combustion  models  which  were  available  prior  to  the  formulation  of  this 
injection  model  required  that  the  user  input  the  fuel-air  ratio  distribution.  This  value  or  values 
was  then  used  in  the  combustion  rate  calculations.  This  requirement  represented  a  restriction 
on  the  validity  of  the  results  in  that  the  fuel  injection  system  response  was  essentially 
assumed. 

The  major  problem  area  with  this  approach  was  augmentor  instabilities  during  transient 
operation  as  the  augmentor  fuel  flowrate  was  increased  towards  maximum  augmentation. 
There  were  numerous  situations  where  rumble  would  occur  on  transients  when  the 
steady-state  fuel  flowrate  was  insufficient  to  cause  rumble.  Sometimes  the  rumble  would 
disappear  as  steady-state  operation  was  achieved  and  other  times  a  rumble-induced  blowout 
would  occur.  This  indicated  that  a  fuel  maldistribution  occurred  which  was  significant  enough 
to  cause  rumble  during  the  transient  operation  of  the  fuel  system. 

The  sprayring  model  was  formulated  to  evaluate  this  transient  operat 'he  app>-  h 
which  was  selected  is  to  specify  the  manner  in  which  fuel  is  introduced  qw  ,  >he  fuel  injector 
system  and  then  evaluate  the  response  of  this  system  to  the  flowrate.  The  mode!  performs 
thermodynamic  and  hydrodynamic  balances  among  the  sprayring  fuel,  the  sprayring  wall 
material,  and  the  specified  control  in  order  to  predict  the  fuel  outflow  distribution  from  the 
injector  at  a  specified  time. 

1'he  type  of  allowed  fuel  control  and  the  specified  fuel  sprayrings  (or  spraybars)  was 
made  as  general  as  possible  to  analyze  a  wide  selection  of  augmentor  system  designs. 

c.  Model  Description 

This  model,  based  on  the  simplified  physical  model  shown  schematically  in  Figure  12,  will 
handle  sprayrings  fed  from  either  one  or  two  locations  or  spraybars  fed  from  one  end.  It  will 
analyze  the  start  transients  in  hot  multiple  zone  systems  with  various  zones  started  simultane¬ 
ously,  sequentially  or  in  some  combination  of  the  two  zones.  The  model  accounts  for  the 
interaction  between  fuel  flow,  heat  transfer,  and  fuel  thermodynamic  characteristics,  predict¬ 
ing  the  instantaneous  distribution  of  thermodynamic  states  throughout  the  fuel  system  and 
the  fuel  flowrate  distribution  into  the  augmentor.  Evaluation  of  an  augmentor  start  transient 
considers  the  fuel  system  initially  hot  and  filled  with  air,  and  analyzes  the  fill  process  as  air 
flows  out  of  the  system  and  the  sprayring  cools  down.  The  model  simultaneously  evaluates  a 
set  of  continuity  relations  for  air,  liquid  and  gaseous  fuel;  a  set  of  energy  relations  for  the  fuel 
and  sprayring,  and  calculates  the  flowrates,  heat  transfer  and  thermodynamic  characteristics 
required  for  solution  of  these  equations. 


18 


(0Aui)  jno  JodeA 


Figure  12.  Schematic  of  Spraybar  Showing  Fluid  Regimes 


The  continuity  relations  are  developed  by  considering  the  mass  rates  of  change  of  air, 
liquid,  and  fuel  vapor  in  the  system,  assuming  all  three  components/phases  are  at  the  same 
pressure,  and  solving  the  resulting  equations  simultaneously  for  the  interdependent  variables. 
The  result  is  a  set  of  rate  equations  defining  the  rates  of  change  of  system  pressure,  liquid 
volume,  and  vapor  volume  in  the  system.  The  resulting  continuity  set,  written  for  a  spraybar 
fed  from  one  end,  is: 


dt  Apv  p„ 


dx  _  da  lx-x„)  m,  -  mlo  -  mv 
dt  dt  pA 

I'lte  dependent  variables  are  system  pressure  (P),  fuel-air  interface  location  (x), 
two  phase  interface  (x„)  and  vapor  void  fraction  (a).  The  relative  flowrates  are  shown  in  Figure 
17.  and  m,  is  the  rate  of  vapor  generation.  Liquid,  vapor  and  air  density  (p,  /\,  pa),  and 
cross-sectional  area  (A),  and  specific  heat  ratio  for  air  (y  are  the  remaining  independent 
variables.  From  the  nature  of  this  equation  set,  it  is  apparent  that  heat  transfer  and  flow 
calculations  must  be  performed  prior  to  or  simultaneously  with  the  evaluation  of  the  set,  and 
are  in  turn  dependent  upon  the  results. 

Figure  FI  shows  the  logic  flow  developed  to  provide  a  satisfactory  approximation  to 
simultaneous  solution  of  the  system  equations  without  the  need  for  costly  iteration. 

Capability  is  incorporated  in  the  model  to  account  for  several  complex  flow  phenomena 
in  the  fuel,  including  two-phase  choking  flow,  flow  with  boiling  dryout,  and  instability  with 
surging  How.  Heat  transfer  models  include  capability  for  analysis  of  heat  transfer  from  the  hot 
gas  stieam  to  the  sprayring  and  heat  transfer  from  the  sprayring  to  the  fuel.  Evaluation  of  heat 
transfer  in  the  fuel  as  the  intially  hot  sprayring  is  cooled  considers  all  possible  modes  of  boiling 
and  convection,  including  film  boiling  with  subcooled  or  saturated  fluids,  nucleate  boiling  with 
subcooled  fuel,  and  convection.  These  models  are  integrated  to  allow  evaluation  of  the 
thermodynamic  change  of  state  from  the  fuel  entry  condition  until  it  is  fully  vaporized. 

The  model  will  analyze  either  single  or  multiple  zone  augmentor  fuel  systems  with  inlet 
condition  specified  as  either  flow  control  at  the  augmentor  entry,  or  inlet  pressure  control  with 
a  variety  of  feedback  parameters.  Transient  analysis  of  the  fuel  system  provides  instantaneous 
calculation  of  system  pressure,  temperature  distributions  in  the  fuel  and  sprayring,  flow 
distributions  into  the  augmentor  and  thermodynamic  states  of  the  fuel  as  it  is  injected,  either 
in  the  vapor  or  liquid  condition. 

Figure  14  shows  the  arrangement  of  a  sprayring  with  fuel  supply  at  two  locations.  Figures 
15  to  17  show  the  results  of  analysis  of  a  start  transient  for  a  sprayring,  with  initial  air-filled 
volume,  for  operation  in  a  flow-controlled  mode.  Figure  15  illustrates  how  system  pressure 
responds  to  the  filling  process  and  Figures  16  and  17  show  the  filling  process  and  fuel 
discharge.  These  results  can  be  used  to  determine  the  fuel  distribution  into  the  augmentor. 
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Figure  13.  Flow  Logic  for  Augmentor  Fuel  System  Model 


Figure  17.  Location  of  Fuel/ Air  Interface  in  Sprayring 


The  basic  sprayring  model  requires  that  the  total  inlet  flow  for  a  sprayring  is  supplied  to 
the  analysis.  The  model  then  generates  the  sprayring  pressure  and  injected  flowrate  distribu¬ 
tion  as  a  function  of  time.  Therefore,  two  options  have  been  written  to  supply  flow  to  the 
system. 

One  such  control  is  a  distribution  valve.  Such  a  valve,  shown  schematically  in  Figure  18, 
is  used  to  sequentially  open  flow  to  each  sprayring  as  a  function  of  time.  The  flow  schedule  is 
shown  in  Figure  19. 
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Figure  IS.  Distribution  Valve  Schematic 


Rather  than  use  a  ramped  flowrate,  it  also  may  be  input  directly.  Either  the  total 
flowrate  or  individual  sprayring  flowrates  may  be  input  as  a  function  of  time.  If  the  total 
flowrate  is  used,  the  percentage  which  goes  into  each  sprayring  must  be  specified. 

The  interface  requirements  between  the  sprayring  model  and  the  flameholder  combustion 
model  are  primarily  those  of  geometric  and  aerodynamic  consistency.  The  output  from  the 
analysis  of  a  given  sprayring  consists  of  fuel  flowrate  distribution  around  the  circumference  of 
the  sprayring.  The  sprayring  is  analytically  divided  into  a  given  number  of  segments,  typically 
20.  Each  segment  has  a  calculated  amount  of  injected  fuel  flowrate  which  occurs  at  the  time 
specified  by  the  user. 

This  fuel  flowrate  output  is  then  divided  among  the  streamtubes  of  the  correct  portion  of 
the  flameholder.  This  relationship  is  shown  in  Figure  20.  The  fuel-air  ratio  in  this  particular 
streamtube  will  be  set  by  the  amount  of  fuel  flowrate  from  the  sprayring  and  the  calculated 
airflow  from  the  sprayring  and  the  calculated  airflow  from  the  flameholder  description.  For 
this  reason,  care  must  be  taken  to  ensure  that  the  input  for  the  sprayring  and  flameholder 
descriptions  are  consistent.  For  example,  the  combustion  model  may  be  successfully  executed 
using  only  one  streamtube  each  to  describe  the  fan  and  core  streams  of  the  augmentor.  If  the 
sprayring  option  is  selected,  however,  all  of  the  fuel  from  a  given  sprayring  would  be  placed 
into  each  single  streamtube,  resulting  in  a  very  rich  fuel-air  ratio. 

The  methodology  by  which  the  correct  assignment  is  made  is  given  in  the  following 
section. 

The  flameholder  description  is  input  as  an  array  of  description  data.  The  first  streamtube 
is  the  one  at  top  dead  center  and  proceeding  clockwise  around  the  flameholder  (looking 
forward  in  the  engine).  Each  sprayring  is  told  the  index  of  the  initial  streamtube  and  the  final 
streamtube  which  it  feeds.  The  program  will  then  calculate  which  section  of  the  sprayring 
feeds  each  streamtube  and  allocate  the  fuel  from  that  section  to  the  given  streamtube.  This  is 
shown  schematically  in  Figure  21. 
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Figure  20.  Sprayring  -  Flameholder  Geometric  Compatibility 

One  final  input  is  required  to  define  the  geometric  compatibility.  The  fuel  flowrates  from 
the  sprayring  are  calculated  referenced  to  the  location  of  the  only  or  primary  feed  location 
while  the  flameholder  is  referenced  to  top  dead  center.  To  finish  the  angular  interface,  the 
amount  of  counterclockwise  offset  of  the  feed  from  top  dead  center  is  input  (see  Figure  211. 
The  fuel-air  ratio  for  each  streamtube  is  calculated  from  this  amount  of  fuel  flowrate  and  the 
airflow  in  the  associated  streamtube.  The  airflow  calculation  is  based  upon  the  inlet  conditions 
to  that  streamtube,  input  in  the  arrays,  and  the  delta  radius  which  the  appropriate  sprayring 
covers.  These  are  shown  as  DRLRU)  and  DELR(2)  in  Figure  21.  If  the  sprayring  and 
combustion  models  are  run  sequentially,  the  arrays  NSC  and  NSH  should  he  set  to  1  to 
maintain  airflow  consistency  in  the  program. 

5.  PROGRAM  SETUP 

The  combined  sprayring/rumble/flameholder  combustion  model  program  supplied  by- 
Pratt  &  Whitney  Aircraft  contains  all  the  subroutines  necessary  to  operate  the  programs  with 
the  exception  of  systems  routines  normally  supplied  by  the  computer  manufacturer.  The 
program  is  written  in  Fortran  IV  and  runs  on  any  large  scale  computer  system  with  little  or  no 
modification  required.  Test  case  input  and  output  are  included  to  verify  successful  installation. 
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6.  PROGRAM  PERFORMANCE  OPTIONS 

The  combined  model  computer  program  has  options  which  control  the  models  which  are 
executed,  as  well  as  the  type  of  sprayring  and  control:  the  augmentor  type,  the  fan  splitter 
type,  the  combustion  data  source,  fuel  type  and  output  format.  These  options  are  explained  in 
the  following  sections: 

a.  Run  Option 

The  program  is  stuctured  to  permit  execution  of  any  one  of  the  three  major  models  or  the 
following  combinations  of  them: 

Input  Symbol  Description 

IRUN  If  IRUN  =  1,  the  program  executes  the  Sprayring  Model 

only. 

If  IRUN  -  2,  the  program  executes  the  Flameholder 
Combustion  Model  only. 

If  IRUN  =  ii,  the  program  executes  the  Rumble  Model 
only. 

If  IRUN  -  4,  the  program  executes  the  Sprayring  Model 
and  Flameholder  Combustion  Model  in  sequence. 

If  IRUN  =  S,  the  program  executes  the  Flameholder 
Combustion  Model  and  the  Rumble  Model  in  sequence. 

If  IRUN  =  6,  the  program  executes  all  three  models  in 
sequence. 

b.  Sprayring  Model  Options 

There  are  two  options  available  which  control  the  description  of  the  augmentor  fuel 
flowrate  which  are  used  in  the  execution  of  the  program.  This  section  describes  the  influence 
of  these  options. 

Control  Option  Description 


ICAL  If  ICAL  =  1,  the  sprayring  model  produces  plots  of 

pressure  vs  time  and  flowrate  vs  time. 

If  ICAL  =  0,  no  plots  are  generated. 

IFLOW  If  IFLOW  =  1,  the  linear  flowrate  ramp  is  used  for  each 

sprayring  as  shown  in  Figure  19. 

If  IFLOW  =  2,  a  single  flowrate  curve  as  a  function  of 
time  is  input  using  NTAB,  TTAB,  FTAB.  The  input 
variable  PSPLIT  can  be  used  as  a  zone  multiplier  in  this 
curve. 
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c.  Rumble  Model  Options 


Augmentor  Option 

The  rumble  model  is  designed  to  simulate  three  augmentor  designs:  V -gutter,  Vorbix  or 
Swirl. 

Input  Symbol  Description 

NAUGOP  If  NAUGOP  =  I,  the  rumble  model  simulates  a  V-gutter 

flameholder  augmentor. 

If  NAUGOP  =  2,  the  rumble  model  simulates  a  Vorbix 
augmentor. 

If  NAUGOP  =  3,  the  rumble  model  simulates  a  Swirl 
augmentor. 


Splitter  Option 

The  rumble  model  is  designed  to  simulate  two  fan  splitter  designs:  proximate  splitter  or 
remote  splitter. 

Input  Symbol  Description 

NFSOP  If  NFSOP  =  1,  the  rumble  model  uses  a  proximate 

splitter  assumption  at  fan  discharge  (Fan  duct  does  not 
communicate  with  core  at  fan  discharge). 

If  NFSOP  =  2,  the  rumble  model  uses  a  remote  splitter 
assumption  at  fan  discharge  (Fan  duct  communicates 
with  core  at  fan  discharge). 

d.  General  Options 
Fuel  Options 

The  combined  model  is  designed  to  operate  with  fuels  of  different  lower  heating  values. 
Input  Symbol  Descriptions 

JFUEL  If  JFUEL  =  1,  the  program  uses  values  for  JP-4  fuel. 

If  JFUEL  =  2,  the  program  uses  values  for  JP-5  fuel. 


Print  Options 


The  rumble  model  provides  either  tabular  and  plotted  output  or  just  plotted  output.  The 
flameholder  combustion  model  provides  either  limited  or  full  tabular  output. 

Input  Symbol  Description 

NPRNTR  If  NPRNTR  =■  0,  the  program  provides  both  tabular 

rumble  model  output  and  Open  Loop  Transfer  Function 
plots. 

If  NPRNTR  =  1,  the  program  provides  only  Open  Loop 
Transfer  Function  plots. 

NPRNTF  If  NPRNTF  =  0,  the  program  provides  limited 

flameholder  combustion  model  tabular  output. 

If  NPRNTF  -  1,  the  program  provides  full  flameholder 
combustion  model  tabular  output. 

NPRINT  The  input  integer  for  NPRINT  controls  the  number  of 

time  steps  between  printouts  of  the  sprayring  model. 


7.  INPUT 

a.  General 

The  combined  model  uses  various  input  parameters  depending  on  which  sprayring, 
combustion  and  augmentor  options  have  been  selected  from  the  previously  described  run 
options.  An  input  data  flow  schematic  is  presented  in  Figure  22.  The  specitic  input  blocks, 
which  are  required  for  each  major  run  option  selected,  are  shown  in  the  chart  in  Figure  23. 
Figures  24  and  25  are  schematics  of  the  rumble  model  and  the  flameholder  combustion  model 
geometry  identification.  The  user  should  note  that  not  all  input  parameters  are  required  for 
any  given  run  option. 


Figure  22.  Combined  Model  Input  Requirements  Attachment 


Figure  24.  Rumble  Model  Geometry  Identification 
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Figure  25.  Flameholder  Combustion  Model  Geometry  Schematic 


b.  Input  Description 

The  program  uses  “Namelist”  input  as  defined  in  the  applicable  computer  manual,  i.e., 
for  the  IBM  370,  the  manual  is  IBM  Systems  360/370  Fortran  IV  Language  Manual, 
GC‘28-6515-10,  pages  54-56.  This  section  gives  the  “Namelist”  names  and  a  description  for 
each  input  variable. 


Parameter  Name 

Description 

ALM(16) 

Sprayring  pintle  area  parameter 

ALPHAC(IOO) 

Fan  stream  flameholder  apex  angle,  degrees 

ALPHAHt  100) 

Core  stream  flameholder  apex  angle,  degrees 

ALPHAO 

Maximum  void  fraction  for  100%  nucleate  boiling  in 
sprayring  (must  be  <1.00) 

AM(  16) 

Sprayring  cross-sectional  flow  area,  sq.  ft. 

AO(16) 

Maximum  flowrate  ramped  flow  if  IFL0W  =  1 

BLMU6) 

Pintle  flow  area  parameter,  dimensionless 
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Parameter  Name 


BPR  Bypass  ratio,  fan  duct  airflow/core  air-flow,  dimensionless 

CF  Film  coefficient  multiplier  for  sprayring  heat  transfer, 

dimensionless 

DEL  Time  increment  for  sprayring  model  calculations,  seconds 

DELR(16)  Delta  radius  of  augmentor  into  which  fuel  flow  from  each 

sprayring  flows,  inches.  Refer  to  Figure  21 

DF  Film  coefficient  Reynolds  number  for  sprayring  heat 

transfer,  dimensionless 

DM(16)  Sprayring  passage  perimeter,  feet 

DPCS  Fan  side  vane  pressure  loss  (AP/P)  from  mixing  plane  to 

ignition  plane  (STA  3  to  STA4),  dimensionless  (Swirl 
augmentor  only) 

DPD  Fan  duct  pressure  loss  (AP/P)  allocated  to  STA  2, 

dimensionless.  Allocate  remainder  to  STA  3;  see  DPS 

DPH  Pressure  loss  (AP/P)  from  mixing  plane  to  ignition  plane 

(STA  3  to  STA  4),  dimensionless.  For  V-gutter  augmen¬ 
tor  this  accounts  for  spraybar  and  flameholder  pressure 
loss.  For  Vorbix  augmentor  this  accounts  for  Vortex  gen¬ 
erator  and  pilot  pressure  loss  (core  and  fan  combined) 

DPHS  Core  side  vane  pressure  loss  (AP/P)  from  mixing  plane  to 

ignition  plane  (STA  3H  to  STA  4),  dimensionless  (Swirl 
augmentor  only) 

DPS  Fan  duct  pressure  loss  (AP/P)  allocated  to  STA  3. 

dimensionless.  Allocate  remainder  to  STA  2;  see  DPD 

DTA(16)  Delta  time  for  control  valve  for  each  sprayring  to  fully 

open,  sec. 

DTC  Maximum  temperature  difference  for  film  boiling  calcula¬ 

tion  in  sprayring,  °R.  (DTC  =  Tw  Ts|) 

DTSTR  Maximum  temperature  difference  for  nucleate  boiling  in 

sprayring,  °R.  (DTSTR  =  Tw  —  Ts|) 


DTO 


Minimum  temperature  difference  for  nucleate  boiling,  °R 


Parameter  Name 
EPSC 
EPSH 
ETA 

ETAC 

ETAH 

FA 


FAC 

FAC(IOO) 

FAH 

FAH(IOO) 

FAV 

FHWC(IOO) 

FHWH(IOO) 

FTAB(50) 

LA 


LB 


LC 


Description 

Fan  stream  turbulence  factor,  dimensionless 

Core  stream  turbulence  factor,  dimensionless 

Augmentor  overall  combustion  efficiency,  actual 
temperature  rise/ideal  temperature  rise,  dimensionless 

Augmentor  fan  stream  combustion  efficiency,  actual 
temperature  rise/ideal  temperature  rise,  dimensionless 

Augmentor  core  stream  combustion  efficiency,  actual 
temperature  rise/ideal  temperature  rise,  dimensionless 

Augmentor  overall  fuel-air  ratio,  dimensionless.  Defined 
as  augmentor  total  fuel  flow/fan  stream  air  flow  (STA  3) 
plus  core  stream  air  flow  (STA  3H)  plus  primary  engine 
fuel  flow  (STA  3H) 

Augmentor  fan  stream  fuel-air  ratio,  dimensionless.  De¬ 
fined  as  augmentor  fan  stream  fuel  flow/fan  stream  air 
flow  (STA  3) 

Augmentor  fuel-air  ratio  for  each  individual  fan  stream 
flow  tube,  dimensionless 

Augmentor  core  stream  fuel-air  ratio,  dimensionless.  De¬ 
fined  as  augmentor  core  stream  fuel  flow/core  stream  air 
flow  (STA  3H)  plus  primary  engine  fuel  flow  (STA  3H) 


Augmentor  fuel-air  ratio  for  each  individual  core  stream 
flow  tube,  dimensionless 

Vitiated  fuel-air  ratio  of  core  stream  at  entry  to  augmen¬ 
tor  (STA  3H),  dimensionless.  Defined  as  primary  engine 
fuel  flow  (STA  3H)/core  stream  air  flow  (STA  3H) 

Individual  flameholder  widths  in  fan  stream,  inches 

Individual  flameholder  widths  in  core  stream,  inches 

Tabular  flowrate  values,  lbm/hr 

Length  of  augmentor,  mixing  plane  to  nozzle  (STA  3  to 
STA  II),  inches 

Distance  from  ignition  plane  to  nozzle  (STA  4  to  STA 
11),  inches 

Length  of  fan  duct,  fan  discharge  to  mixing  plane  (STA  1 
to  STA  3),  inches 
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Parameter  Name 


Description 


LH  Distance  from  turbine  discharge  to  mixing  plane  (STA 

2H  to  STA  3H),  inches 

LI  Distance  from  ignition  plane  to  beginning  of  combustion 

zone  (STA  4  to  STA  5),  inches 

LK  Distance  from  ignition  plane  to  end  of  combustion  zone 

(STA  4  to  STA  10),  inches 


LM(16)  Sprayring  heated  length,  ft. 

LSC  Distance  from  spraybar  to  flameholder  in  fan  stream, 

inches 


LSC(IOO)  Distance  from  spraybar  to  flameholder  for  each  individual 

streamtube  in  the  fan  stream,  inches 

LSH  Distance  from  spraybar  to  flameholder  in  core  stream, 

inches 


LSH(IOO)  Distance  from  spraybar  to  flameholder  for  each  individual 

streamtube  in  the  core  stream,  inches 

L2  Distance  from  fan  duct  pressure  loss  (DPD)  to  mixing 

plane  (STA  2  to  STA  3),  inches 

M6C  Fan  stream  Mach  number  at  entry  to  augmentor  (STA  3), 

dimensionless.  (Must  be  >0.) 


M6H  Core  stream  Mach  number  at  entry  to  augmentor  (STA 

3H),  dimensionless.  (Must  be  0.) 

M6R  Mach  number  of  mixed  augmentor  stream  flow  prior  to 

combustion  (STA  4),  dimensionless.  (Must  be  0.) 

NOST(16)  Number  of  streamtubes  fed  by  each  sprayring,  integer 

(Must  be  >1.) 


NSC(IOO) 

NSH(IOO) 

NSR 

NSRC 

NSRH 

NTAB 


Number  of  fan  stream  flow  tubes  of  this  type,  integer 

Number  of  core  stream  flow  tubes  of  this  type,  integer 

Total  number  of  sprayrings,  integer 

Total  number  of  sprayrings  feeding  fan  air  flow,  integer 

Number  of  sprayrings  which  feed  the  core  stream  air  flow, 
integer 

Number  of  points  in  flowrate  table 

Number  of  streamtube  types  in  the  fan  flow,  integer 
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NTC 


Parameter  Name 
NTH 

OFFSET*  16) 
PCTM(16) 

PFSR(IOO) 

PO 

PRNOZ 

PSPLIT(16) 

PS6 

SBLF(16) 

STBTIM 

TA(16) 

TAUCOM 

TAUC(IOO) 

TAUH(IOO) 

TCORE 

TEXT*  100) 


Description 

Number  of  streamtube  types  in  the  core  flow,  integer 

Angle  of  sprayring  feed  line  relative  to  top  dead  center  of 
augmentor,  degrees.  Refer  to  Figure  21  for  details 

Length  split  between  4  and  4  directions  of  sprayring 
relative  to  feed  line,  dimensionless  (must  be  <1.0) 

Individual  spraybar  fuel  pressure  for  each  fan  flow  tube, 
psia 

Initial  sprayring  internal  pressure  prior  to  fuel  flow.  Gen¬ 
erally  equal  to  PS6,  psia 

Exhaust  nozzle  pressure  ratio  (always  >1.),  dimension¬ 
less.  If  nozzle  is  choked,  any  value  greater  than  critical 
value  required  to  choke  nozzle  (approximately  2.0)  may 
be  input.  Exact  value  of  PRNOZ  is  required  only  if  nozzle 
is  unchoked 

Fraction  of  input  total  flowrate  fed  to  a  given  sprayring 
(must  be  1.0),  dimensionless 

Augmentor  static  pressure  at  entry  to  augmentor  (STA 
3),  psia 

Length  of  a  spraybar  which  is  exposed  to  fan  airflow,  feet 

Time  at  which  stability  calculation  is  performed,  sec 

Time  at  which  the  control  valve  for  a  given  sprayring 
starts  to  open,  seconds 

Average  value  of  flameholder  wake  residence  time,  sec¬ 
onds 

Individual  streamtube  blockage  ratio  for  fan  stream, 
dimensionless 

Individual  streamtube  blockage  ratio  for  core  stream, 
dimensionless 

Core  engine  time  constant,  mass  of  air  in  core  engine 
volume/mass  flowrate  of  air  through  the  core  engine,  sec 

External  flow  temperature  for  individual  flow  tubes  in  the 
fan  flow,  °R 
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Parameter  Name 


TF 

TFHMT(IOO) 

TFSR(IOO) 

THM<16) 

TIME1 

TIME2 

TMAX 

TO 

TOL 

TTAB(50) 

TWC(16) 

TWF(16) 

T3H 

T6C(100) 

T6H(100) 

V1M(16) 

WCOOL 

WEXTOOO) 


1 


Description 

Fuel  temperature  at  inlet  to  control,  °R 

Flameholder  metal  temperature  for  each  streamtube  in 
fan  stream,  °R 

Sprayring  fuel  temperature  for  individual  streamtubes  in 
fan  stream,  °R 

Sprayring  wall  thickness,  feet 

First  time  for  stability  analysis  of  combined  models, 
seconds 

Second  time  for  stability  analysis  of  combined  models, 
seconds 

Maximum  time  for  sprayring  calculations,  seconds 

Initial  air  temperature  within  sprayring,  °R 

Iteration  between  tolerance  in  sprayring  model, 
dimensionless 

Times  for  flowrate  table,  sec 

Initial  wall  temperature  of  portion  of  spraybar  in  core 
stream,  °R 

Initial  wall  temperature  of  portion  of  spraybar  in  fan 
stream,  °R 

Main  burner  inlet  temperature,  °R 

Fan  stream  temperature  at  entry  to  augmentor  (STA  3), 
°R 

Core  stream  temperature  at  entry  to  augmentor  (STA 
3H),  for  individual  flow  tubes,  °R 

Volume  of  sprayring  ahead  of  split  in  feed  lines,  cubic  ft 

Ratio  of  nozzle  cooling  air  to  total  engine  air  flow, 
dimensionless 

External  flow  ratio  for  individual  flow  tubes  in  the  fan 
stream,  dimen'  -nless 
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Parameter  Name 


XLC(IOO) 

XLH(IOO) 

ZEF 


ZEFC 


ZEFH 


ZEFP 


ZEP 


ZEPC 


ZEPH 


Description 

Distance  from  flameholder  to  nozzle  for  individual  fan 
stream  flow  tubes,  inches 

Distance  from  flameholder  to  nozzle  for  individual  core 
stream  flow  tubes,  inches 


Normalized  slope,  augmentor  overall  combustion  effi¬ 
ciency  vs  overall  fuel-air  ratio, 


FA  (tETA 
ETA  UFA 


dimensionless 


Normalized  slope,  augmentor  fan  stream  combustion  effi¬ 
ciency  vs  fan  stream  fuel-air  ratio, 

FAC  (1ETAC 

ETAC  ~1faT  ’  dimensl"nless 


Normalized  slope,  augmentor  core  stream  combustion  ef¬ 
ficiency  vs  core  stream  fuel-air  ratio, 


FAH  oETAH 
ETAH  <9FAH 


dimensionless 


Normalized  slope,  augmentor  overall  combustion  effi¬ 
ciency  vs  fuel-air  ratio  of  the  pilot  burner, 


FAP  dETA 
ETA  dFAP 


dimensionless 


Normalized  slope,  augmentor  overall  combustion  effi¬ 
ciency  vs  pressure  at  ignition  plane, 


P  r/ETA 
ETA  ^P 


dimensionless 


Normalized  slope,  augmentor  fan  stream  combustion  effi¬ 
ciency  vs  pressure  at  ignition  plane, 


P  dETAC 
ETAC  it  P 


dimensionless 


Normalized  slope,  augmentor  core  stream  combustion  ef¬ 
ficiency  vs  pressure  at  ignition  plane, 


P  dETAH 
ETAH  dP 


dimensionless 
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Parameter  Name 

ZETC 

ZETH 

ZEVC 

ZEVH 


Description 

Normalized  slope,  augmentor  fan  stream  combustion  effi¬ 
ciency  vs  fan  stream  entry  temperature. 


T6C  dETAC 


ETAC  ST6C 


,  dimensionless 


Normalized  slope,  augmentor  core  stream  combustion  ef¬ 
ficiency  vs  core  stream  entry  temperature, 


T6H  dETAH 


ETAH  dT6H 


,  dimensionless 


Normalized  slope,  augmentor  fan  stream  combustion  effi¬ 
ciency  vs  fan  stream  entry  velocity, 

V  dETAC  .  , 

- -  ,  dimensionless 

Normalized  slope,  augmentor  core  stream  combustion  ef¬ 
ficiency  vs  core  stream  entry  velocity, 


V  dETAH  ..  .  . 

etah'  — dV —  ’  “,mensionless 
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c.  Input  Setup 


In  addition  to  the  “Namelist”  input,  the  program  requires  input  for:  (a)  additional 
optional  ratio  calculations,  (b)  output  plot  selection  and  format  ,  and  (c)  frequency  range  and 
increment  selection.  The  input  setup  is  shown  in  Figure  150  and  is  described  below. 


(1)  Each  input  case  requires  a  title  card.  Column  1  for  the  first  case  must 
contain  a  1.  The  “Namelist”  input  must  be  preceded  by  an  &  INPUT 
starting  in  column  2  and  followed  by  an  &  END  starting  in  column  2.  The 
“Namelist”  input  (Columns  2-80)  required  for  each  case  is  presented  in 
Figures  22  and  23.  Each  input  must  be  separated  by  a  comma  (see  Figure 
26).  The  ratio  calculations,  plot  setup  and  frequency  selection  must  follow 
the  first  input  case.  For  additional  input  cases,  follow  the  frequency 
selection  cards  with  a  blank  card  and  then  the  additional  title  cards  and 
input  cases.  Only  those  parameters  that  differ  from  the  previous  case 
must  be  input.  For  the  additional  input  cases,  if  column  1  of  the  title  card 
contains  a  1,  the  ratio  calculations,  plot  setup  and  frequency  selection  will 
be  the  same  as  the  preceding  case.  If  column  1  of  the  title  card  contains 
a  1,  new  ratio  calculations,  plot  setup  and  frequency  selection  may  be 
input. 

(2)  Additional  ratio  calculations  may  be  performed  by  inputting  the  parame¬ 
ter  identification  numbers  as  indicated  in  Figure  26.  Up  to  40  ratios  may 
be  calculated  and  these  ratios  will  automatically  be  included  in  the 
tabular  output.  The  parameter  identification  numbers  are  presented  in 
Figure  27.  One  blank  field  will  terminate  this  type  input.  If  columns 
71-75  are  used,  a  blank  card  must  follow. 


(3)  Calcomp  plots  of  any  rumble  model  output  parameter  may  be  obtained 
by  inputting  one  card  for  each  parameter  as  described  below.  A  maximum 
of  10  plots  may  be  requested  for  any  case.  A  blank  card  must  be  input  to 
terminate  plot  requests  or  if  no  plots  are  desired. 


Column  1  -  3 

Column  11  -  20 
Column  21  -  30 
Column  31  -  40 
Column  41  -  50 
Column  51-60 
Column  61  -  70 
Column  71-80 


—  Output  Parameter  No.,  right  adjusted 
(integer;  no  decimal) 

—  Amplitude  Option  (decimal  required) 

—  Phase  Option  (decimal  required) 

—  Frequency  Option  (decimal  required) 

—  Amplitude  Factor  (decimal  required) 

—  Frequency  Factor  (decimal  required) 

—  XMIN  (decimal  required) 

—  XMAX  (decimal  required) 
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FORTRAN  CODING  FORM 

_ _  ENGINEER 


figure  26.  Combined  Rumble  Model  In  pul  Setup 


Output  Parameter  Parameter  Identification 
Name  Number 


Pt 

1 

VI 

2 

R1 

3 

P2 

4 

V2 

5 

R2 

6 

P3 

7 

V3 

8 

R3 

9 

P3H 

10 

V3H 

1 1 

R3H 

12 

P2H 

13 

V2H 

14 

R2H 

15 

QIN 

16 

W3 

17 

W3H 

18 

QOUT 

19 

P4 

20 

V4 

21 

R4 

22 

P5 

23 

V5 

24 

R5 

25 

P6 

26 

V6 

27 

R6 

28 

P7 

29 

V7 

30 

R7 

31 

P8 

32 

V0 

33 

R8 

34 

P9 

35 

V9 

36 

R9 

37 

P10 

J8 

V10 

39 

RIO 

40 

P11 

41 

VII 

42 

R11 

4? 

Note:  PI  =  P1/Qin  =  (APi/Pi)/(AQ1N/Q)N)  Same  tor  Other  Output  Parameters 


Figure  27.  Parameter  Identification  Numbers 
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Output  Parameter  No. 


• 

A  list  of  output  parameter  identification  numbers  is  presented  in  Figure 
31.  If  any  of  the  additional  ratio  calculations  (described  above)  are  to  be 
plotted,  parameter  identification  numbers  starting  at  101  and  in¬ 
cremented  by  1  are  used.  The  parameter  identification  numbers  must  be 
right  adjusted  in  columns  1-3. 

•  Amplitude  Option 

If  Amplitude  Option  =  0.,  Log  (amplitude)  will  be  plotted. 

If  Amplitude  Option  =  1.,  Amplitude  will  be  plotted. 

•  Phase  Option 

If  Phase  Option  =  ().,  Phase  angle  of  0  to  -360  will  be  plotted. 

If  Phase  Option  =  1.,  Phase  angle  of  180  to  -180  will  be  plotted. 

•  Frequency  Option 

If  Frequency  Option  —  0.,  Log  (Frequency)  will  be  plotted. 

If  Frequency  Option  =  I.,  Frequency  will  be  plotted. 

•  Amplitude  Factor 

If  Amplitude  Option  =■  0.,  the  Amplitude  Factor  must  be  input  as  a 
number  which  could  be  written  as  10N,  where  N  =  ( ± )  integer.  Log 
(Amplitude  Factor)  will  be  added  to  the  base  Log  amplitude  scale,  where 
the  base  Log  amplitude  scale  ranges  from  -1.0  to  2.0. 

If  Amplitude  Option  =  L,  the  Amplitude  Factor  becomes  a  multiplier  for 
the  base  amplitude  scale,  where  the  base  amplitude  scale  ranges  from  0. 
to  3.0. 

•  Frequency  Factor 

If  Frequency  Option  =  0.,  the  Frequency  Factor  must  be  input  as  a 
number  which  could  be  written  as  10N,  where  N  =  ( ± )  integer.  Fre¬ 
quency  factor  will  be  a  multiplier  for  the  base  frequency  scale,  where  the 
base  frequency  scale  ranges  from  .1  to  100. 

If  Frequency  Option  =  I.,  the  Frequency  Factor  will  be  a  multiplier  for 
the  base  frequency  scale,  where  the  base  frequency  scale  ranges  from  0.  to 
100.,  unless  XMIN  and  XMAX  are  input.  XMIN  is  the  minimum  value 
for  the  frequency  scale  when  Frequency  Option  =  1.  XMAX  is  the 
maximum  value  for  the  frequency  scale  when  Frequency  Option  =  1.  If 
either  Amplitude  Factor  or  Frequency  Factor  is  input  at  0.  when  XMIN 
and  XMAX  are  input,  the  program  automatically  sets  Amplitude  Factor 
or  Frequency  Factor  to  1. 
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(4)  The  frequencies  used  in  the  program  calculations  are  input  in  two  parts. 

First,  there  are  three  cards  which  contain  the  minimum  frequency, 
increment  and  maximum  frequency  (see  Figure  24).  The  increment  is 
used  to  determine  each  frequency  for  the  range  defined.  Additional 
independent  frequencies  (up  to  500  values)  may  be  input  in  fields  of  10, 

7  fields  per  card.  A  field  containing  -1.  will  terminate  this  input.  If  the 
7th  field  of  the  last  card  is  used,  an  additional  card  with  -1.  in  the  first 
field  is  required. 

(5)  The  rumble  model  has  been  set  up  to  model  a  turbofan  engine  with  a 
mixed  flow  augmentor.  To  model  a  turbojet  (no  fan),  set  BPR  =  0  and 
NFSOP  =  1.  To  model  a  fan  duct  augmentor  (separate  fan  and  core 
flows),  set  BPR  =  1010  and  NFSOP  =  1. 

8.  OUTPUT 

a.  General 

This  section  is  presented  in  three  parts:  (1)  sprayring  model  output,  (2)  combustion 
model  output,  and  (3)  rumble  model  output.  The  tabular  or  printed  output  for  all  output 
options  is  presented  for  the  three  models.  The  outputs  are  presented  in  the  order  in  which 
they  appear  for  a  full-combined  model  run  (IRUN  =  6). 

b.  Output  Description 

(1)  The  first  output  is  a  full  printing  of  the  input  for  the  case  which  is  run. 

This  is  titled  INPUT  LISTING  and  is  a  direct  repeat  of  all  input  values. 

(2)  The  second  output  is  a  printing  of  all  default  parameters.  If  a  required 
input  value  has  not  been  listed,  the  program  will  supply  a  reasonable 
value  from  an  internal  library.  The  user  should  check  this  listing  to 
determine  if  the  default  value  is  appropriate  to  the  particular  problem 
being  analyzed. 

(3)  The  third  section  prints  all  the  namelist  parameters  in  an  expanded 
format  for  verification  purposes.  This  list  includes  both  input  and  default 
values  and  represents  the  actual  parameter  values  used  in  the  computer 
run. 

(4)  The  next  section  will  write  any  check  errors  for  erroneous  input  values.  If 
an  input  error  is  detected,  the  program  will  terminate. 
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(5)  If  the  sprayring  and  flameholder  combustion  model  have  both  been 
executed  (IRUN  =  4  or  6),  the  next  output  section  will  print  the  results 
of  the  geometric  compatibility  subroutine.  This  subroutine  (ASSOC) 
assigns  a  specific  portion  of  one  sprayring  to  a  given  streamtube.  Some  of 
the  output  is  an  appropriate  input  listing.  The  new  output  parameters 
are: 


Parameter  Name 

SR 

Description 

Sprayring  identification  number 

NOST 

Number  of  streamtubes  fed  by  this  spray¬ 
ring 

NFST 

The  number  of  the  first  streamtube  fed  by 
this  sprayring 

NLST 

The  number  of  the  last  streamtube  fed  by 
this  sprayring 

SUMDEG 

Total  degrees  of  sprayring  to  be  analyzed. 
Equal  to  360  and  OFFSET 

SECDEG 

Number  of  degrees  in  the  repeating  spray- 
ring  sectors 

TUBDEG 

Number  of  degrees  of  sprayring  per  stream- 
tube 

ST 

Streamtube  identification  number 

IBEGH 

IBEGC 

Identification  number  of  the  first  sprayring 
sector  feeding  this  streamtube  in  fan 
(IBEGC)  or  core  (IBEGH). 

IENDH 

IENDC 

Identification  number  of  the  last  sprayring 
sector  feeding  this  streamtube  in  fan 
(IENDC)  or  core  (IENDH). 

FRACBH 

FRACBC 

The  fraction  of  IBEGC/IBEGH  which  feeds 
this  streamtube 

FRACEH 

FRACEC 

The  fraction  of  IENDC/IENDH  which  feeds 
this  streamtube 

WIDTH 

The  width  of  this  streamtube,  inches. 
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(6)  Sprayring  Model  Output 

This  section  presents  the  results  of  the  sprayring  model.  A  review  of  the 
appropriate  input  is  printed  first.  Next  follows  a  tabular  listing  of  the 
model  output  as  a  function  of  time.  The  time  between  print  steps  is 
controlled  by  DEL  and  NPRINT  inputs. 


Parameter  Name 
TIME 


Description 

Time  from  the  start  of  the  flow  calculations, 
seconds 


N 

A(N) 

P(l) 

X 


X-XO 


Y 

Y-YO 

ALPX 

ALPY 

PRESS 

FLOWIN 

TW(I) 

TW(E) 


Identification  number  of  the  sprayring 
Distribution  valve  open  area,  ft* 

Supply  pressure,  psia 

Axial  location  of  the  fluid  fuel  face  in  the 
x-direction,  feet 

Length  of  the  two-phase  region  in  the 
x-direction,  feet 

Axial  location  of  the  fluid  fuel  face  in  the 
y-direction,  feet 

Length  of  the  two-phase  region  in  the 
y-direction,  feet 

Void  fraction  in  the  x  two-phase  region 
Void  fraction  in  the  y  two-phase  region 
Sprayring  pressure,  psia 
Fuel  flowrate  into  the  sprayring,  lbm/hr 
Wall  temperature  at  sprayring  inlet,  °R 
Wall  temperature  at  sprayring  exit,  °R 


TOTFLOW  Total  fuel  flow  rate  out  of  the  sprayring  into 

the  streamtubes,  lbm/hr 

L  FLOW  Liquid  fuel  flowrate  out  of  the  sprayring 

into  the  streamtubes,  lbm/hr 


AXMXO  Equals  ALPX  •  (X-XO) 


AYMXO  Equals  ALPY  (Y-YO) 

At  the  completion  of  the  tabular  output  versus  time,  the  calculated  fuel 
flowrate  and  fuel -air  ratio  values  are  printed.  The  fuel  flowrate  is  printed  in 
Ibm/sec  units.  Also  printed  are  the  fuel  flowrates  from  each  section  of  the 
sprayring  (SRFUEL)  in  Ibm/sec  units. 
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(7)  Flameholder  combustion  model  full  tabular  output. 


Parameter(s) 

Fan  Stream 

Streamtube  type 
No.  of  this  type 

Static  Pressure  (PS6) 

Approach  Temperature  (T6C) 
Approach  Mach  No.  (M6C) 
Input  FA  Ratio  (FAC) 

Effective  FA  Ratio 

F/H  Width  (FHWC) 

Blockage  Ratio  (TAUC) 

F/H  Apex  Angle  (ALPHAC) 
S/R  Fuel  Temperature  (TFSR) 


Description 

Identifies  following  sections  as  fan  duct  out¬ 
put 

Identifies  for  this  set  of  input  variables 

The  number  of  streamtubes  with  this  set  of 
input  variables 

Inlet  static  pressure,  psia 

Inlet  temperature,  °R 

Inlet  flow  Mach  No.,  d’less 

Inlet  fuel-air  ratio,  d’less 

Effective  fuel-air  ratio  accounting  for  liner 
cooling  air  flow 

Flameholder  width,  inches 

Ratio  of  flameholder  width  to  stream  tube 
width,  d’less 

V-gutter  flameholder  apex  angle,  degrees 

Temperature  of  the  fuel  within  sprayring, 
°R 


S/R  Fuel  Pressure  (PFSR)  Pressure  of  the  fuel  within  the  sprayring, 

psia 

S/R  to  F/H  Distance  (LSC)  Axial  separation  distance  between  the 

sprayring  and  the  flameholder,  inches 

F/H  to  Nozzle  Distance  (XLC)  Axial  distance  from  the  flameholder  to  the 

exhaust  nozzle  throat,  inches 


Turbulence  Level  (EPSC)  Ratio  of  RMS  turbulence  velocity  to  the 

approach  velocity  at  the  inlet,  d’less 

Wake  Flow  Addition  (WEXT)  Ratio  of  external  wake  flow  to  recirculated 

flow,  d’less 


Flow  Source  Temperature  Temperature  of  above  flow,  °R 

(TEXT) 

Effective  Inlet  Temperature  Mass  average  of  WEXT  flow  at  TEXT  and 

recirculated  flow  at  T6C,  °R 
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Parameters 


Description 


Chemical  combustion  efficiency  Average  efficiency  baser!  on  average  stream- 

tube  exit  temperature  and  average  effective 
fuel-air  ratio,  d’less 

Thermal  combustion  efficiency  Average  efficiency  based  on  streamtube  av¬ 
erage  exit  temperature  plus  cooling  air  and 
average  input  fuel-air  ratio,  d’less 


Average  cooling  air 
temperature 


Mass  averaged  inlet  temperature  used  for 
cooling,  °R 


Average  streamline  exit 
temperature 


Mass  average  of  the  streamtubes  without 
cooling  air,  °R 


Average  duct  exit  temperature  Mass  average  of  streamtubes  plus  cooling 

air,  °R 


Total  flowrate 


Total  of  each  streamtube  type  times  the 
number  of  each  type,  Ibm/sec 


Average  fuel  air-ratio 
Core  stream 


Mass  average  of  the  input  fuel-air  ratios 

Identified  following  sections  as  core  stream 
output 


Wake  recirculation  coefficient  Same  as  Kl  in  fan  duct,  d’less 


Ideal  temperature  rise 


M/S  Fuel-air  ratio 


Ideal  temperature  rise  based  on  input 
fuel-air  ratio  and  main  burner  fuel-air  ratio. 
See  Appendix  B 

Fuel-air  ratio  of  the  vitiated  air  entering  the 
core  streamtubes 


M/B  Inlet  temperature 


Inlet  temperature  to  the  main  burner,  d’less 


Average  distance  from  spraybar  Average  axial  distance  from  the  spraybars  to 
to  F/H  the  flameholders,  inches 

Note:  Any  core  stream  parameters  which  are  not  listed  above  have  the  same 
definition  as  their  fan  stream  counterpart. 

(4)  Flamehoider  combustion  model  limited  tabular  output. 


Parameters 
Fan  stream 


Description 

Identifies  following  as  fan  stream  cases  in¬ 
put  and  output 


Inlet  temperature  (T6C) 


Streamtube  inlet  temperature,  °R 


Fuel-air  ratio  (FAC) 


Input  fuel-air  ratio,  d’less. 
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Parameter(s) 


Description 


Average  inlet  temperature 


Mass  averaged  inlet  temperature,  °R 


Ideal  temperature  rise 


Streamtube  ideal  temperature  rise  based  on 
effective  fuel-air  ratio,  °R 


Combustion  efficiency 
Exit  temperature 


Streamtube  combustion  efficiency,  d’less 

Streamtube  exit  temperature  based  on  effec¬ 
tive  fuel-air  ratio,  °R 


Average  ideal  temperature  rise  Ideal  temperature  rise  based  on  average  in¬ 
put  fuel-air  ratio,  °R 

Average  combustion  efficiency  Efficiency  based  on  average  ideal 

temperature  rise  with  cooling  air  effect  in¬ 
cluded,  d’less 


Average  exit  temperature 


Core  stream 


Exit  temperature  including  cooling  air,  °R 
Identifies  following  as  core  stream  section 


Mach  No. 


Streamtube  inlet  flow  Mach  number,  d’less 


Average  ideal  temperature  rise  Ideal  temperature  rise  based  on  average  in¬ 
put  fuel-air  ratio,  °R 

Average  combustion  efficiency  Efficiency  based  on  average  temperature  rise 

and  average  ideal  temperature  rise,  d’less 


Average  exit  temperature 


Exit  temperature  based  on  mass  averaged 
actual  temperature  rise,  °R 


(8)  Rumble  model  tabular  output  (listed  in  the  order  they  appear): 


Parameter(s) 


Description 


NAMELIST  INPUT 


KNOZ 


The  “namelist”  input  parameters  and  the 
values  input  are  listed  for  verification 

A  parameter  that  relates  the  influence  of 
pressure  at  STA  11  on  velocity  at  ST  A  11, 
dimensionless 


FAAB 


Augmentor  overall  fuel-air  ratio,  dimension¬ 
less 


ETAAB 


DTIAB 


Augmentor  overall  efficiency,  dimensionless 

Augmentor  overall  ideal  temperature  rise, 
°R 
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Parameter(s) 


Description 


DTAB 


Augmentor  overall  actual  temperature  rise, 
°R 


T6M 


Augmentor  mixed  temperature  before  com¬ 
bustion  (ST A  3),  °R 


TKC 

XLHV 

DTC 

QCQT 

DTIC 

TAUDC 

DTH 

QHQT 

DTIH 

TAUDH 

ZTFC 


ZTFH 


L  (1-11) 
YL  (1-11) 


Augmentor  mixed  exhaust  temperature 
(STA  10),  °R 

Lower  heating  value  for  the  fuel  selected, 
Btu/lbm 

Fan  stream  temperature  rise,  °R 

Fraction  of  total  heat  release  contributed  to 
fan  stream,  dimensionless 

Fan  stream  ideal  temperature  rise,  °R 

Fan  stream  drift  delay  from  spraybar  to 
flameholder,  sec 

Core  stream  temperature  rise,  °R 

Fraction  of  total  heat  release  contributed  bv 
core  stream,  dimensionless 

Core  stream  ideal  temperature  rise,  °R 

Core  stream  drift  delay  from  spraybar  to 
flameholder,  sec 


Normalized  slope,  augmentor  fan  stream 
ideal  temperature  rise  vs  fan  stream  fuel-air 
ratio, 

FAC  ./DTIC  ,. 

-  -  ,  dimensionless. 

DTIC  dFAC 


DTIC  FAC 


Normalized  slope,  augmentor  core  stream 
ideal  temperature  rise  vs  core  stream  fuel-air 
ratio 

FAH  3 DTIH 


dimensionless. 


DTIH  <9FAH 


Distance  between  model  stations,  in. 

Station  locations  references  to  STA  1,  inches 
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Parameter(s) 

Description 

C  (1-11) 

Velocity  of  sound  at  each  station,  in. /sec 

CH 

Velocity  of  sound  in  core  stream  at  STA  3H, 
in./sec 

M  (1-11) 

Mach  number  at  each  station,  dimensions 

MH 

Mach  number  in  core  stream  at  STA  3H, 
dimensionless 

T  (1-11) 

Temperature  at  each  station,  °R 

TH 

Temperature  in  core  stream  at  STA  3H,  °R 

PRHOT 

Pressure  drop  through  combustion  zone 
(STA  5  -  STA  10),  psia 

G  (1-11) 

Ratio  of  specific  heats  at  each  station, 
dimensionless 

GH 

Ratio  of  specific  heats  in  a  core  stream  at 
STA  3H,  dimensionless 

TAUF  (1-11) 

Time  delays  for  downstream  running  sonic 
waves  between  stations,  sec 

TAUFH 

Time  delay  for  downstream  running  sonic 
wave  between  STA  2H  and  3H,  sec 

TAUG  (1-11) 

Time  delays  for  upstream  running  sonic 
waves  between  stations,  sec 

TAUGH 

Time  delay  for  upstream  running  sonic  wave 
between  STA  2H  and  3H,  sec 

TAUE  (1-11) 

Time  delays  for  downstream  running  en¬ 
tropy  waves  between  stations,  sec 

TAUEH 

Time  delay  for  downstream  running  entropy 
wave  between  STA  2H  and  3H,  sec 

QOP  (1-11) 

Ratio  of  volumetric  heat  release  rate  at  each 
station  to  pressure  at  each  station,  1/sec 

(2)  Rumble  model  Open  Loop  Transfer  Function  plots.  Each  plot  consists  of 
Open  Loop  Transfer  Function  Gain  vs  frequency  and  phase  vs  frequency 

for  the  parameters  selected. 
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9.  PROGRAM  MESSAGES  AND  LIMITS 


a.  Input  C hacks 

The  program  checks  all  inputs  to  see  if  the  inputs  are  missing  or  equal  the  default  values 
built  into  the  deck.  Missing  inputs  are  set  equal  to  the  default  values.  A  warning  message 
(presented  below)  is  printed  to  alert  the  user  if  default  input  values  are  identified.  The 
program  also  checks  specific  inputs  to  ensure  reasonable  input  data.  If  these  data  checks  are 
not  satisfied,  the  run  will  be  canceled.  These  checks  and  corresponding  print-out  messages  are 
presented  below: 


Condition 

Input  value  =  default  value  or  no 
input  for  certain  parameter 

LA  <  LB,  where  LK  and  LB  are 
Rumble  Model  Inputs 

LK  >  LB,  where  LK  and  LB  are 
Rumble  Model  Inputs 

LA  <  LCALC,  where  LCALC  =  LB 
+  MAX  (LSC,  LSH) 

LI  >  LK 

LC  <  L2 

ETA  <  0  or  >  1. 

ETAC  <  0  or  >  1. 

ETAH  <  0  or  >  1. 

(FAC  +  FAH)  =  0 


Message 

WARNING  —  PARAMETER 
XXXXXX  =  YYY.Y  is  a  default 
value 

INPUT  ERROR  1  —  LA  must  be 
greater  than  or  equal  to  LB 

INPUT  ERROR  2  —  LB  must  be 
greater  than  or  equal  to  LK 

INPUT  ERROR  3  -  LA  must  be 
greater  than  or  equal  to  the  sum  of 
LB  plus  the  max  of  LSC  or  LSH.  LA 
has  been  adjusted  accordingly. 
Check  input. 

INPUT  ERROR  4  —  LI  must  be 
less  than  LK. 

INPUT  ERROR  5  -  LC  must  be 
greater  than  or  equal  to  L2. 

INPUT  ERROR  6  —  ETA  must  be 
between  0  and  1. 

INPUT  ERROR  7  —  ETAC  must 
be  between  0  and  1. 

INPUT  ERROR  8  —  ETAH  must 
be  between  0  and  1. 

INPUT  ERROR  9  —  FAC  and  FAH 
cannot  both  be  zero  with  augmentor 
on. 
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C 


H 


Condition 

[FAV  +  (1  +  FAV)  FAH1 

r  XLHVi 
l  18500.  J  ~ 

II?An>  r  XLHV  1  ^ 
or  if  (FAC)  18500  J  _ 

NFSOP  -  2  and  BPR  =  0 

DPCS  <  0  or  >  1. 

DPD  <  0  or  >  1. 

DPH  <  0  or  >  1. 

DPHS  <  0  or  >  1. 

DPS  <  0  or  >  1. 

T3H  >  2200. 

BPR  <  0. 

FAV  <  0. 

NAUGOP  <  0.  or  >  3. 

NFSOP  <  0  or  >  2. 

IRUN  <  0  or  >  6. 

JFUEL  <  0  or  >  2. 


Men. sat;*' 

x  INPUT  ERROR  10  —  Core  or  fan 

stream  totai  fuel-air  ratio  exceeds 
limits  of  ideal  temperature  rise 
curve.  Blowout  likely. 


0.09 


INPUT  ERROR  1 1  -  BPR  cannot 
be  zero  when  the  remote  flow  split¬ 
ter  option  is  selected. 

INPUT  ERROR  12.  -  DPCS  must 
be  between  0  and  1. 

INPUT  ERROR  13  —  DPD  must  be 
between  0  and  1. 

INPUT  ERROR  14  -  DPH  must  be 
between  0  and  1. 

INPUT  ERROR  15  -  DPHS  must 
be  between  0  and  1. 

INPUT  ERROR  16  -  DPS  must  be 
between  0  and  1. 

INPUT  ERROR  17  —  T3H  exceeds 
limits  of  ideal  temperature  rise 
curve.  T3H  must  be  less  than 
2200°  R. 

INPUT  ERROR  18  —  BPR  must  be 
equal  to  or  greater  than  0. 

INPUT  ERROR  19  —  FAV  must  be 
equal  to  or  greater  than  0. 

INPUT  ERROR  20  —  NAUGOP 
must  be  1,  2  or  3. 

INPUT  ERROR  21  —  NFSOP  must 
be  1  or  2. 

INPUT  ERROR  22  —  IRUN  must 
be  1,  2,  3.  4,  5  or  6. 

INPUT  ERROR  23  —  JFUEL  must 
be  1  or  2. 
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Condition 

Message 

NPRNTR  <  0  or  >  1. 

INPUT  ERROR  24  -  NPRNTR 
must  be  0  or  1. 

NPRNTF  <  0  or  >  1. 

INPUT  ERROR  25  —  NPRNTF 
must  be  0  or  1. 

M6C  <  0 

INPUT  ERROR  26  —  M6C  must  be 
greater  than  0. 

M6H  <  0 

INPUT  ERROR  27  —  M6H  must 
be  greater  than  0. 

M6R  <  0 

INPUT  ERROR  28  —  M6R  must  be 
greater  than  0. 

LI  <  0 

INPUT  ERROR  29  —  LI  must  be 
equal  to  or  greater  than  0. 

LK  <  0 

INPUT  ERROR  30  —  LK  must  be 
greater  than  0. 

LA  <  0 

INPUT  ERROR  31  -  LA  must  be 
greater  than  0. 

LB  <  0 

INPUT  ERROR  32  -  LB  must  be 
greater  than  0. 

LC  <  0 

INPUT  ERROR  33  -  LC  must  be 
greater  than  0. 

LSC  >  LA 

INPUT  ERROR  34  —  LSC  must  be 
less  than  LA. 

LSH  >  LA 

INPUT  ERROR  35  -  LSH  must  be 
less  than  LA. 

LH  <  0 

INPUT  ERROR  36  —  LH  must  be 
greater  than  0. 

L2  <  0 

INPUT  ERROR  37  -  L2  must  be 
greater  than  or  equal  to  0. 

TCORE  <  0 

INPUT  ERROR  38  —  TCORE  must 
be  equal  to  or  greater  than  0. 

PS6  <  0 

INPUT  ERROR  39  —  PS6  must  be 
greater  than  0. 

T6C  <  0 

INPUT  ERROR  40  —  T6C  must  be 

greater  than  0. 
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Condition 


Messant 


T6H  <  0 

T3H  <  460. 

FA  <  0 

FAC  <  0 

FAH  <  0 

PRNOZ  <  1 

ALPHAC  (100)  < 

ALPHAH  (100)  < 

ESPC  <  0  or  >  1 

EPSH  <  0  or  >  1 

FAC  (100)  <  0 

FAH  (100)  <  0 

FAV  >  0.068 

FHWC  (100)  <0 
FHWH  (100)  <  0 


INPUT  ERROR  41  -  T6H  must  be 
greater  than  0. 

INPUT  ERROR  42  —  T3H  must  be 
greater  than  460. 

INPUT  ERROR  43  FA  must  be 
greater  than  or  equal  to  0. 

INPUT  ERROR  44  —  FAC  must  be 
greater  than  or  equal  to  0. 

INPUT  ERROR  45  —  FAH  must  be 
greater  than  or  equal  to  0. 

INPUT  ERROR  46  —  PRNOZ  must 
be  greater  than  1. 

0  or  >  180  INPUT  ERROR  47  ALPHAC 

must  be  greater  than  0  and  less  than 
or  equal  to  180  deg. 

0  or  >  180  INPUT  ERROR  48  —  ALPHAH 

must  be  greater  than  0  and  less  than 
or  equal  to  180  deg. 

INPUT  ERROR  49  —  EPSC  must 
be  greater  than  or  equal  to  0  and 
less  than  or  equal  to  1. 

INPUT  ERROR  SO  —  EPSH  must 
be  greater  than  or  equal  to  0  and 
less  than  or  equal  to  1. 


INPUT  ERROR  51  —  FAC  must  be 
greater  than  0. 

INPUT  ERROR  52  -  FAH  must  be 
greater  than  0. 

INPUT  ERROR  53  —  FAV  cannot 
exceed  stoichiometric  (0.068). 

INPUT  ERROR  54  —  FHWC  must 
be  greater  than  0. 

INPUT  ERROR  55  —  FHWH  must 
be  greater  than  0. 

INPUT  ERROR  56  —  ESC  must  be 
greater  than  0. 


LSC  (100)  <  0 
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Condition 


Message 


LSH  (100)  <  0 

M6C 

1-TAUC  (100) 

M6H 


>  1. 


1-TAUH  (100)'  >  L 
NSC  (100)  <  0  or  >  100. 


NSH  (100)  <  0  or  >  100. 


NTC  <  0  or  >  100. 


NTH  <  0  or  >  100 


•  *00)  PS6 

TAUC  (100)  <0or>  1. 

TAUH  (100)  <0or>  1. 

TEXT  (100)  <  460. 

TFSR  (1(H))  <  460. 

T6C  (100)  <  460. 

T6H  (100)  <  460. 

WEXT  (100)  <0 


INPUT  ERROR  57  —  LSH  must  be 
greater  than  0. 

INPUT  ERROR  58  —  Flow  is  su¬ 
personic  in  fan  stream  at  the 
flameholder  plane. 

INPUT  ERROR  59  Flow  is  su¬ 
personic  in  core  stream  at  the 
flameholder  plane. 

INPUT  ERROR  60  —  NSC  must  be 
greater  than  or  equal  to  0  and  less 
than  or  equal  to  100. 

INPUT  ERROR  61  —  NSH  must  be 
greater  than  or  equal  to  0  and  less 
than  or  equal  to  100. 

INPUT  ERROR  62  -  NTC  must  be 
greater  than  or  equal  to  zero  and 
less  than  or  equal  to  100. 

INPUT  ERROR  63  -  NTH  must 
be  greater  than  or  equal  to  0  and 
less  than  or  equal  to  100. 

INPUT  ERROR  64  -  PFSR  must 
be  greater  than  PS6. 

INPUT  ERROR  65  —  TAUC  must 
be  greater  than  0  and  less  than  1. 

INPUT  ERROR  66  —  TAUH  must 
be  greater  than  0  and  less  than  1. 

INPUT  ERROR  67  —  TEXT  must 
be  greater  than  or  equal  to  460. 

INPUT  ERROR  68  —  TFSR  must 
be  greater  than  or  equal  to  460°R. 

INPUT  ERROR  69  —  T6C  must  be 
greater  than  460. 

INPUT  ERROR  70  —  T6H  must  be 
greater  than  460. 

INPUT  ERROR  71  -  WEXT  must 
be  greater  than  or  equal  to  0. 
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Condition 
XLC  (100)  <  0 

XLH  (100)  <  0 

WCOOL  <  0.  or  >  1. 

ALM  (16)  <  0 

ALPHAO  <  0  or  >  1 

AM  (16)  <  0 
AO  (16)  <  0 

BLM  (16)  <  0 

CF  <  0 

DEL  <  0 

DELR  (16)  <0 

DF  <  0 

DM  (16)  <  0 

If  IFLOW  =  1,  DTA  (16)  <  0 


Message 

INPUT  ERROR  72  —  XLC  must  be 
greater  than  0. 

INPUT  ERROR  73  —  XLH  must 
be  greater  than  0. 

INPUT  ERROR  74  —  WCOOL 
must  be  greater  than  or  equal  to  0 
and  less  than  1. 

INPUT  ERROR  75  —  ALM  must 
be  greater  than  0. 

INPUT  ERROR  76  —  ALPHO 
must  be  greater  than  or  equal  to  0 
and  less  than  1. 

INPUT  ERROR  77  —  AM  must  be 
greater  than  0. 

INPUT  ERROR  78  —  AO  must  be 
greater  than  0  if  IFLOW  =  1. 

INPUT  ERROR  79  —  BLM  must 
be  greater  than  or  equal  to  0. 

INPUT  ERROR  80  —  CF  must  be 
greater  than  0. 

INPUT  ERROR  81  —  DEL  must  be 
greater  than  0. 

INPUT  ERROR  82  —  DELR  must 
be  greater  than  0. 

INPUT  ERROR  83  —  DF  must  be 
greater  than  or  equal  to  0. 

INPUT  ERROR  84  —  DM  must  be 
greater  than  0. 

INPUT  ERROR  87  —  DTA  must  be 
greater  than  or  equal  to  0,  if 
IFLOW  =  1. 
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Condition 

If  IFLOW  =  1,  DTA  (16)  >  TMAX 

DTC  <  0 

DTSTR  <  0 

DTO  <  0 

DTO  >  DTSTR 

DTSTR  >  DTC 

LM  (16)  <0 

NOST  (16)  <  1 

NSR  <  1 

NSRC  <  0 

NSRH  <  0 

NSR  =  (NSRC  +  NSRH) 

OFFSET  (16)  <  0  or  >  180 

PCTM  (16)  >  1.0 


Message 

INPUT  ERROR  88  —  DTA  must  be 
less  than  TMAX,  if  IFLOW  =  1. 

INPUT  ERROR  89  -  DTC  must  be 
greater  than  0. 

INPUT  ERROR  90  —  DTSTR  must 
be  greater  than  0. 

INPUT  ERROR  91  —  DTO  must  be 
greater  than  0. 

INPUT  ERROR  92  —  DTO  must  be 
less  than  DTSTR. 

INPUT  ERROR  93  —  DTSTR  must 
be  less  than  DTC. 

INPUT  ERROR  95  —  LM  must  be 
greater  than  0. 

INPUT  ERROR  96  —  NOST  must 
be  greater  than  or  equal  to  1. 

INPUT  ERROR  97  —  NSR  must  be 
greater  than  or  equal  to  1. 

INPUT  ERROR  98  -  NSRC  must 
be  greater  than  or  equal  to  0. 

INPUT  ERROR  99  —  NSRH  must 
be  greater  than  or  equal  to  0. 

INPUT  ERROR  100  -  NSR  must 
equal  NSRC  plus  NSRH. 

INPUT  ERROR  101  —  OFFSET 
must  be  greater  than  or  equal  to  0 
and  less  than  or  equal  to  180. 

INPUT  ERROR  102  —  PCTM  must 
be  less  than  or  equal  to  1. 
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Condition 

Message 

PO  <  0 

INPUT  ERROR  105  —  PO  must  be 
greater  than  0. 

PSPLIT  (16)  <0 

INPUT  ERROR  106  —  PSPLIT 
must  be  greater  than  or  equal  to  0. 

SBLF  (16)  <  0 

INPUT  ERROR  109  —  SBLF  must 
be  greater  than  0. 

STBTIM  <0or>  TMAX 

INPUT  ERROR  110  —  STBTIM 
must  be  greater  than  0  and  less  than 
or  equal  to  TMAX. 

TA  (16)  <  0  or  >  TMAX 

INPUT  ERROR  111  —  TA  must  be 
greater  than  or  equal  to  0  and  less 
than  or  equal  to  TMAX. 

TAUCOM  <  0 

INPUT  ERROR  112  —  TAUCOM 
must  be  greater  than  or  equal  to  0. 

TF  <  0 

INPUT  ERROR  113  —  TF  must  be 
greater  than  0. 

TFSR  (100)  <0 

INPUT  ERROR  114  —  TFSR  must 
be  greater  than  0. 

THM  (16)  <  0 

INPUT  ERROR  115  —  THM  must 
be  greater  than  0. 

TMAX  <  0 

INPUT  ERROR  116  —  TMAX 
must  be  greater  than  0. 

TO  <  0 

INPUT  ERROR  117  —  TO  must  be 
greater  than  0. 

TOL  <  0  or  >  TMAX 

INPUT  ERROR  118  —  TOL  must 
be  greater  than  or  equal  to  0  and 
less  than  TMAX. 

TWC  (16)  <  0 

INPUT  ERROR  119  —  TWC  must 
be  greater  than  0. 

TWF  (16)  <  0 


INPUT  ERROR  120  —  TWF  must 
be  greater  than  0. 


Condition 
VIM  (16) 

V2M  (16) 
V3M  (16) 

If  IFLOW 

IFLOW  * 


1 


Message 


0 

INPUT  ERROR  121  —  VlM  must 
be  greater  than  0. 

0 

INPUT  ERROR  122  —  V2M  must 
be  greater  than  0. 

0 

INPUT  ERROR  123  —  V3M  must 
be  greater  than  0. 

2,  NTAB  <  0  or  >  50 

INPUT  ERROR  129  —  NTAB  must 
be  greater  than  0  and  less  than  or 
equal  to  50  when  IFLOW  =  2. 

or  2 

INPUT  ERROR  132  —  IFLOW 
must  be  equal  to  0  or  1. 
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Condition  Message 


If  IFLOW  =  2,  PSPLIT  (16)  <  0  INPUT  ERROR  136  —  PSPL1T 

must  be  greater  than  0  when 
IFLOW  =  2. 


NPRINT  <  1  INPUT  ERROR  138  —  NPRINT 

must  be  equal  to  or  greater  than  1. 

NPRINT  >  (TM AX/DEL)  INPUT  ERROR  139  —  NPRINT 

cannot  exceed  the  number  of  time 
steps  calculated  (TMAX/DEL). 


b.  Calculation  Failure  Messages  for  Flameholder  Combustion  Model 

The  Flameholder  Combustion  Model  program  checks  for  specific  calculation  failures.  The 
causes  and  corresponding  messages  are  presented  below: 


Message 

Warning  ***  wake  temperature  iter¬ 
ation  failed  for  streamtube  No.  XX 


Aerodynamic  loading  exceeds  kinetic 
capacity 


All  fuel  vaporized  —  terminate  case 


Overall  fuel-air  ratio  below  lean 
limit 


FAR  outside  flammability  limit 


Cause 

The  calculated  value  of  the  fan  duct 
flameholder  wake  fuel-air  ratio  ex¬ 
ceeds  the  rich  limit  at  the  inlet  con¬ 
ditions  for  this  streamtube.  This 
streamtube  case  has  failed. 

The  inlet  values  of  velocity,  pres¬ 
sure,  and  temperature  produce  a 
wake  loading  which  exceeds  the  re¬ 
action  limit  at  this  fuel-air  ratio. 
This  streamtube  case  has  failed. 

The  injection  process  has  resulted  in 
only  vapor  fuel.  Run  this  input  as  a 
core  case  if  desired. 

The  input  fuel-air  ratio  is  less  than 
the  calculatd  minimum  value  for 
flame  propagation  in  the  fan  stream. 

The  input  fuel-air  ratio  is  outside 
the  limits  of  the  data  for  laminar 
flame  speed  built  into  the  program. 
Currently  set  at  0.027  lean  and  0.120 
rich  limit  fuel-air  ratios. 


10.  PROGRAM  LISTINGS 

Complete  listings  of  the  computer  program  are  provided  in  Appendix  C. 
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11.  TEST  CASES 


r\ 


The  test  cases  which  follow  are  provided  in  Appendix  D.  These  cases  were  selected  to 
provide  the  prospective  user  with  a  broad  range  of  the  potential  combinations  of  the  various 
input  options.  All  input  parameters  have  been  used,  even  if  the  value  was  the  same  as  a  default 
value.  Prior  to  executing  a  new  case,  the  user  should  become  familiar  with  the  appropriate  test 
case. 


Test 

Models  Executed 

Case 

Spravrine 

Combustion 

Rumble 

l 

X 

— 

— 

2 

— 

X 

3 

— 

— 

X 

4 

— 

— 

X 

5 

— 

— 

X 

6 

X 

X 

— 

7 

— 

X 

X 

8 

X 

X 

X 

I  RUN 

N  AUG  OP 

Notes 

1 

— 

2 

— 

3 

1 

V -gutter  F/H 

3 

2 

Vorbix  augmentor 

3 

3 

Swirl  augmentor 

4 

— 

5 

1 

6 

1 

All  of  these  were  executed  with  JP4  fuel,  input  fuel  flowrates  for  the  sprayring  model  and 
remote  splitter  for  the  rumble  model.  Full  output  was  selected  in  all  cases. 
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12.  PROGRAM  IDENTIFICATION  AND  REVISION  PROCEDURE 
a.  CCD  Number 

Customer  Computer  Decks  (CCD’s)  are  identified  *•  ^  CCD  number  and  date.  An 

example  is  CCD  1001=0.0  November  15,  1969. 

The  CCD  number  consists  of: 

•  Basic  number  (first  four  digits) 

•  Dash  (no  change)  number 

•  Decimal  (or  addition/correction)  number. 

(1)  Basic  Number 

This  four-digit  number  generally  corresponds  to  a  given  program.  If  another  method  is 
developed  or  studied  which  does  not  replace  or  supersede  the  original,  a  new  four -digit  number 
is  issued. 

(2)  Dash  Number 

Major  changes  in  the  program  are  reflected  in  different  dash  numbers.  A  change  in 
techniques  or  mathematical  methods  would  produce  a  new  dash  number  provided  the  new 
techniques  replace  or  change  the  old  techniques.  For  more  than  nine  changes,  the  dash  number 
shifts  to  letters.  The  dash  number  of  the  original  program  is  zero. 

(3)  Decimal  Number 

The  decimal  is  used  for  all  other  program  changes  each  as: 

•  Adding  new  optional  routines 

•  A  change  to  the  FORTRAN  source  language  due  to  computer  differences 

•  Correcting  a  mistake  in  the  program 

•  Input  and  output  changes 

•  Adding  unique  customer  oriented  curves  or  initialization  data. 

In  most  cases  a  decimal  number  change  is  made  and  documented  by  writing  an 
addendum  to  the  user’s  manual  without  reprinting  the  user’s  manual.  An  errata  or  addendum 
to  the  user’s  manual  is  used  in  conjunction  with  a  Computer  Simulation  Change  Notice  (Figure 
28).  This  notice  briefly  describes  the  changes  affecting  the  deck  and  manual. 

Accompanying  the  notice  is  a  new  manual  title  page  reflecting  the  new  date  and/or  dash 
or  decimal  changes.  The  SCN  also  includes  the  change  pages  with  black  bars  in  the  right-hand 
margin  opposite  data  changed  and  a  new  date  in  the  upper  right-hand  corner. 

Programs  revised  by  the  user  without  written  consent  of  the  supplier  are  the  responsi¬ 
bility  of  the  user. 
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Engine  or  Component  Designation 


Date 


User’s  Manual  No.  _  Customer  Computer  Deck  No.  _ 

Deck  and  User's  Designation  is  Changed  to: _ _ — 

Dated:  _ 

Change  Nature:  |  |  Errata  f  1  Addendum  Changes  Effect:  1  |  Deck  |  )  Manual 

Change:  Reason: 


Detailed  Reasons  for  Changes  are  Shown  In  Enclosure 

Approval: 

Program  Office  _ 

Systems  Stability  and  Control  - 

C.  H.  Borgmeyer 


fO 


Figure  28.  Computer  Simulation  Change  Notice 
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APPENDIX  A 

DEVELOPMENT  TUMBLE  MODEL  EQUATIONS 


1.  DEVELOPMENT  OF  ACOUSTIt.  ATIONS 


In  this  section  equations  are  developed  to  describe  how  velocity,  pressure  and  density 
at  every  point  in  the  augmentor  respond  to  a  combustion  disturbance,  which  is  treated  as  a 
heat  input  to  a  flowing  invicid  ideal  gas  stream.  Knowing  how  these  three  parameters 
(velocity,  pressure,  density)  respond  allows  calculation  of  any  other  parameter  needed,  such 
as  mass  flowrate  or  temperature.  The  first  equations  to  be  developed  are  the  three  longitu¬ 
dinal  wave  equations,  which  are  applicable  between  boundaries  and  discontinuities.  Then 
equations  for  the  boundaries  and  discontinuities  are  developed.  The  wave  equations  plus 
the  boundary  and  discontinuity  equations  are  referred  to  as  the  “acoustic”  equations.  The 
“combustion”  equations  needed  to  complete  the  rumble  model  are  developed  in  paragraph 
2  of  this  appendix. 


Symbols  used  below  are  defined  in  the  list  of  symbols.  For  any  section  of  augmentor  with 
rigid  walls  and  constant  cross-sectional  area,  such  as  shown  in  Figure  29,  through  which  an 
invicid  fluid  (viscosity  is  zero)  is  flowing,  the  one-dimensional  momentum,  continuity  and 
energy  equations  are: 


*  ,pv^  + 
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aV  |  y  AP 
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ax 


,  au 

+  "  « 


(1) 


For  an  ideal  gas,  these  equations  reduce  to  the  following  non-linear  wave  equations: 
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The  wave  equations  are  linearized  by  the  small  perturbation  substitutions: 


P(x.t) 

p(x.t) 

C(x,t) 

V(x.t) 

q(x,t) 


P<  x  I  +  AP(x,t) 
p(x)  +  Ap(x,t) 
C(x)  +  AC(x,t) 
V(x)  +  AV(x,t) 
q(x)  +  Aq(x.t) 


(3) 


Second  order  terms  are  neglected  in  making  the  substitutions. 
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End  Temp 


Fan  Fan  Duct 
Discharge  AP/P  Loss 


Fan  Side 
Mixing 

Plan* 


Fan  Duct 


Core  Engine 


Temp  Rise 


V  -  gutter,  Swirl  or  — '  i  i 
Vorbix  Hardware  (2h)(3h) 


-Combustion  - 
Zone 


Core  Side 
Mixing  Plane 

Turbine 

Discharge 

Figure  29.  Rumble  Model  Station  Identification 


FO  1464*2 

7631 10 
M2- 760 


To  simplify  notation,  the  following  substitutions  are  made  which  normalize  the  change 
n  each  variable  by  its  steady  -state  value: 

i>  =  JJL  v  =—  ^  a  =  (i > 

P  ’  V  •  P  p  '  q  q 

i  In'  ltm-ari/.od  version  of  equations  (2)  becomes: 

+  P+>^V  +  P’  +  >MV'j  +  (7-I)  -p-dr  =  (7- 1) -p  Q' 

P  ym]  +  P  >MV]  +  ~  Cr-U-pq-  if)) 


*i[p  +i[p'-^]  ’  -p*:  =(T-u  J 


3f  =  --  "-2)  [p’  (1-M-M2)  +  PM  +  V  1  +  }  M-M2  [1  +  (1+M)  -|-J  J  J 


Sc,  =  - r p’  (1-M-M2)  +  />M  +  V  j  -i-  +  4-  M-M2  [1  +  (1+M)  -2-1  (  1 

-M' 


/"I 


Taking  the  Laplace  transform  with  respect  to  time,  with  zero  initial  conditions,  and 
letting  subscripts  1  and  2  stand  for  the  upstream  and  downstream  stations  respectively 
(see  Figure  33),  the  general  solution  to  equations  (5)  becomes: 


J 


dx 


8  J 


dx 


„  V+C  A  i  V+C 

(P,  +  -yM.V.I  e  -  (Pi+'yM.V,)  +  .1  dF(x.s)  —  e 


dx 


(>-i)  r  q  „  %  d 
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*  dx 

8  l  v+c 

e  dx 
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(y-  1 


c  q  (l 

s  J  p  A;,x-S)  e 


O  v-c 
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dx 


dx 
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(>-1)  c  q  ,  d 

--j-  /  p  <"«>  -JT 


v-c 

e  dx 


p  dx 

'It 


IP.  -  7P.I  e  -  IP1-7P.I  +  /  4 


J“ 

J  v 


dE(*.s)  —  e 


dx 


(>-l)  f  q  ..  ,  d 

‘  —  1  P  dT  ' 


x  dx 

‘1  V 


dx 


In  equations  (7)  the  first  equation  describes  downstream  running  sonic  waves  of  the  form 
P'  +  V’,  traveling  at  sonic  speed  plus  through-flow  velocity.  The  second  equation  de¬ 

scribes  upstream  running  sonic  waves  of  the  form  P'  -  7M  V’.  traveling  at  sonic  speed  minus 
through-flow  velocity.  The  third  equation  describes  entropy  waves,  p'  -  yp  \  drifting  down¬ 
stream  at  through-flow  velocity. 

The  entropy  waves  become  more  apparent  from  the  expression  for  the  entropy  of  an 
ideal  gas: 

^  =  S  =  P  -  yp'  (8) 

The  entropy  waves  are  related  to  temperature  by: 

yT  =  S'  +  (7-I)  P 

It  is  through  equation  (9)  that  the  drifting  hot  and  cold  combustion  products,  or  entropy 
waves,  are  accounted  for  in  the  rumble  model.  Temperature  changes  produced  as  the 
entropy  waves  strike  the  exhaust  nozzle  create  waves  which  then  travel  back  upstream  at 
sonic  speed. 
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Equations^  7 )  are  not_useful  until  the  integrals  are  evaluated,  which  will  require  defini¬ 
tions  of  V(x),  C(X),  q(x),  P(x)  and  some  assumptions  that  will  allow  integration  of  q  (x.  s), 
dp  (x,  s),  0q  (x,  s)  and  lip  (x,  s).  To  complete  the  solution,  the  augmentor  is  divided  into 
several  “short”  sections,  each  of  length  V,  for  each  of  which  it  can  be  assumed: 


(a) 

d  P(x) 
dx 

0 

tb) 

d  T(x) 
dx 

constant 

/  '• 

to 

q  (x,t) 

-  q  [  0.  t  -  J 

<d) 


qtx) 

P(x) 


constant 


The  small  static  pressure  drop  in  an  augmentor  justifies  assumption  (a).  A  linear  temperature 
rise  throughout  a  section  of  length  1  is  a  good  approximation,  which  justifies  assumption  (b). 
Assumption  (c)  is  the  equation  for  a  “drifting  burning  particle”  releasing  heat  at  a  constant 
volumetric  rate  as  it  drifts  down  the  augmentor.  A  more  detailed  explanation  of  this 
assumption  will  be  provided  in  part  2  (Development  of  Combustion  Equations).  To  justify  the 
constant  steady-state  heat  release  rate  (q),  consider  the  steady-state  version  of  the  energy 
equation  (third  in  equations  (2)). 


V 


1_  dP 
P  dx 


i  s!e. 

J)  dx  - 


With  appropriate  substitutions,  the  equation  reduces  to: 

1  /_>_')  V  dP 

P  '  *,  1  '  P  A  dx  P  dx 


Since 


dP 

dx 


-  0  and 


dT 

dx 


constant,  then 


q_ 

P 


constant 
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.) 


(1U| 


For  a  “short”  section  of  length  V,  the  integration  of  (ip  (x,  s)  in  equations  (7  I  can  In-  carried 
out  as  follows: 


X  _ 

/  -p-  dplx-S)  ^  e  dx  *  ^  0F(o,S)  f  -jjj-  e  dx 
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dx 
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Similar  treatment  allows  integration  of  Pq  (x,  s)  and  p£  (x,  s)  in  equation  17).  To  determine 
how  “short”  a  section  must  be  for  the  solution  to  be  valid,  the  resulting  rumble  model  was 
exercised  repeatedly  while  decreasing  the  section  length  (by  adding  more  stations  in  the 
combustion  zone).  As  the  section  length  decreases,  the  result  will  rapidly  approach  an  exact 
solution.  It  was  found  that  section  lengths  shorter  than  about  20  inches  were  unnecessary. 

With  the  above  assumptions,  equation  (7)  becomes: 
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For  convenience  in  programming  equations  (11)  on  the  computer,  the  following  identity 
substitutions  were  made: 


,  0\{  =  PP,  P  ,  +  RF.p  ,  +  VF,V, 

0  ,2  =  PFj  P  2  +  RFaP  2  +  VF2V2 
0',;,  =  FO,  P ,  +  RGlP  ,  +  VG,V, 
0-„,  =  PG2  P  ,  +■  RGjPl  +  VG2V/ 


(11) 
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whore  by  definition: 
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The  time  constants  in  equations  (12)  were  evaluated  based  upon  the  steady-state 
through-flow  and  sonic  speed  profiles  created  by  the  linear  temperature  gradient. 
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Then  the  time  constants  in  equations  (12)  become: 
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This  completes  the  development  of  the  wave  equations. 

Equations  (11 )  are  applied  throughout  the  augmentor  between  any  two  stations 
between  which  there  is  no  discontinuity.  In  applying  the  equations,  the  general  subscripts  1 
and  2  are  replaced  by  the  actual  upstream  and  downstream  station  numbers,  respectively. 
Referring  to  Figure  33,  they  are  applied  between  stations  (1)  •  (2),  (2)  -  (3),  (4)  -  (5), 

(5)  -  (10)  and  (10)  -  (11).  Between  stations  (1)  through  (5)  and  between  stations  (10)  -  (11 ) 
there  is  no  heat  addition,  and  so  the  heat  addition  terms  q/P  are  set  to  zero.  The  heat  ad¬ 
dition  terms  for  the  combustion  zone,  stations  (5)  -  (10),  are  discussed  in  paragraph  2  of  this 
appendix. 

Discontinuities  occur  at  the  pressure  drop  locations,  stations  (2)  and  (3).  These  are 
modeled  as  small  incompressible  resistive  pressure  drops  of  zero  length.  The  continuity  and 
energy  equations  are  also  applied.  Referring  to  Figure  33,  across  a  pressure  drop: 


P,  -  P,  - 

W,  =  W, 
T,  -  T, 


P7  V  2 


W 


(18) 
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The  equations  are  linearized  and  normalized  as  before  to  yield: 


P, 


U+2Vi) 


pi  +  Vi  -  p',  +  V; 
P«  ~  Pj  =  Ps  —  P» 


(19) 


In  applying  equations  (19)  to  a  given  pressure  drop,  the  general  subscripts  2  and  3  are  re¬ 
placed  by  the  actual  upstream  and  downstream  station  numbers,  respectively.  For  con¬ 
venience  in  programming,  equations  (19)  were  combined  with  the  wave  equations  (1 1 )  to 
eliminate  the  need  for  two  stations  at  each  pressure  drop.  It  is  the  combined  equations 
which  appear  in  the  rumble  model  listing. 


A  junction  occurs  where  the  core  stream  and  fan  stream  enter  the  augmentor  and  form 
the  overall  augmentor  stream  (stations  (3),  (3H)  and  (4)  ).  Again  applying  continuity, 
momentum  and  energy : 

W,  +  W,h  =  W, 
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FAN  SIDE 
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(20) 


W,T,  +  W.„  T,„  =  W,T. 

The  linearized  and  normalized  versions  become : 
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(21) 
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For  the  Swirl  augmentor,  the  momentum  equations  at  stations  (3)  -  (4)  and  (3H)  4) 

are  modified  to  account  for  the  possibility  of  different  pressure  drops  across  the  fan  and  con- 
swirl  vanes.  The  linearized  version  of  the  momentum  equations  for  the  Swirl  augmentor 
becomes: 


Definition  of  the  upstream  and  downstream  boundary  conditions,  at  the  fan  and  at  the  nozzle, 
respectively,  will  complete  the  acoustic  equations.  The  fan  was  assumed  to  be  delivering  a 
constant  mass  flowTate  through  the  fan  OD  (defined  as  that  portion  of  the  fan  between  the  fan 
splitter  and  fan  tip)  and  through  the  fan  ID  (defined  as  that  portion  of  the  fan  between  the 
centerline  and  the  fan  splitter).  It  was  also  assumed  that  the  temperature  of  the  fan  discharge 
flow  could  be  taken  as  time  invarient  (also,  because  of  the  low  Mach  number  at  fan  discharge, 
total  and  static  temperatures  can  be  used  interchangeably).  To  account  for  the  presence  of  a 
core  engine,  and  explore  any  possible  attendant  interaction  with  fan  duct  acoustics,  a  simple 
first  order  lag  representation  of  the  core  engine  was  incorporated  into  the  rumble  model.  The 
core  engine  was  represented  as  a  compressor  delivering  constant  corrected  air  flow  (corrected 
to  compressor  face  conditions)  into  a  lumped  volume.  Flow  out  of  the  volume  exited  through 
a  choked  turbine  to  emerge  at  station  (3H).  The  resulting  transfer  function  for  the  core  engine 
is: 


W,h  =  1 

Pc  1  +  rc  <, he  S 

When-: 

W  =  mass  flowrate  at  station  (3H) 

P,  -  static  pressure  at  the  compressor  face 

r,  (,he  =  core  engine  time  constant 


( 1 


A  default  value  of  tCOre  =  005  seconds  is  built  into  the  rumble  model.  A  different  value  can 
be  input  by  the  user,  and  is  calculated  as  the  mass  of  air  in  the  core  engine  volume  divided  by 
the  mass  flowrate  of  air  through  the  core  engine.  Proximity  of  the  fan  splitter  to  fan  discharge 
also  affects  the  boundary  condition  at  the  fan.  Two  cases  were  considered  and  are  built  into 
the  rumble  model  (see  NFSOP).  In  the  first  case,  called  the  “proximate”  splitter  configuration, 
the  fan  splitter  is  assumed  to  be  so  close  to  fan  discharge  that  no  communication  can  occur 
between  the  fan  duct  and  the  core  engine  across  the  fan  splitter.  For  this  case,  the  boundary 
condition  at  the  fan  becomes: 


P  ,  =  W  -  0 
w,  =  p\  +  V,'  =  0 
T,  =  P,  -  p,'  =  0 
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In  the  second  case,  called  the  “remote”  splitter  configuration,  the  fan  splitter  is  as¬ 
sumed  to  be  sufficiently  remote  from  fan  discharge  to  allow  perfect  communication  between 
the  fan  duct  and  the  core  engine  across  the  fan  splitter.  For  this  case,  the  boundary  condition 
at  the  fan  becomes: 


Pc  =  P, 

W,  =  p,  +  v,-  =  -  BPR 

T,  =  P,  -  p,  =  0 


W, 


P. 


Pl_ 


1  +  rcoHE  S 


(25) 


This  completes  the  definition  of  the  upstream  boundary  condition.  It  is  of  interest  to 
note  that  entropy  waves  are  created  by  sonic  wave  reflections  at  the  upstream  boundary. 
Since  an  entropy  perturbation  is  S'  =  P'.  -  -yp'  ,  and  at  the  boundary  p',  =  P'  ,  then 
S'.  =  (1  -  7)  P'j .  A  similar  argument  will  show  that  entropy  waves  are  also  created  at  the 
pressure  drops  (stations  (2)  and  (3)  ).  These  are  automatically  accounted  for  in  the  rumble 
model,  but  are  of  minor  importance  compared  to  the  entropy  waves  created  in  the  com¬ 
bustion  zone  by  combustion  disturbances. 


The  downstream  boundary  condition  is  based  upon  the  presence  of  a  “short”  nozzle 
just  downstream  of  station  (11),  for  which: 


Wx/TJR 
A  P„ 


0>  (Pr) 


(26) 


where : 


P  R  2> 


P„  =  P0  /nozzle  throat  static  pressure 


When  linearized,  the  downstream  boundary  condition  becomes: 


V„  =-|-  (Pu'-p,. )  +  (KNOZ)  P„ 


where: 


(27) 
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It  is  also  of  interest  to  note  that  for  choked  flow. 


then  KNOZ  =  0  and: 

v„  =|(P,;  -  Pn)  =yTu  (28) 

substituting  from  equation  (16): 

v“  -27s"  +i^p"  1291 

This  equation  directly  relat?s  how  entropy  waves,  as  well  as  pressure  disturbances,  striking  a 
choked  nozzle  will  produce  a  velocity  disturbance. 

This  completes  the  acoustic  equation  development.  These  equations  describe  the 
response  of  pressure,  velocity  and  density  throughout  the  augmentor  to  a  disturbance  in 
combustion.  Development  of  the  corresponding  combustion  equations,  which  describe  how 
combustion  throughout  the  augmentor  will  respond  to  disturbances  in  pressure,  velocity  and 
density,  is  presented  in  the  following  section. 
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2.  DEVELOPMENT  OF  COMBUSTION  EQUATIONS 

Development  of  the  combustion  equations  for  the  V-gutter  flameholder  augmentor  is 
presented  first.  Then  the  combustion  equations  for  Vorbix  and  Swirl  augmentors  are 
presented. 

For  the  V-gutter  flameholder  augmentor  two  combustion  streams,  the  fan  stream  and 
the  core  stream,  are  treated.  This  is  necessary  to  be  able  to  account  for  the  different  com¬ 
bustion  characteristics  of  the  fan  and  core  streams.  The  two  streams  can  have  different 
flameholder  designs  and  fuel-air  ratios  as  well  as  different  flameholder  approach  temperatures 
and  velocities,  causing  the  two  streams  to  have  different  efficiency  vs.  fuel-air  ratio  charac¬ 
teristics.  In  addition,  the  fan  stream  is  preceded  by  a  long  fan  duct  which  can  exhibit  lon¬ 
gitudinal  resonance  at  the  low  frequencies  associated  with  rumble.  The  core  stream  is  pre¬ 
ceded  by  a  short  section  terminating  at  turbine  discharge,  which  is  much  less  responsive  at 
low  frequencies. 

The  basic  approach  taken  for  the  rumble  model  was  to  model  combustion  disturbances 
in  the  fan  and  core  streams  independently,  accounting  for  the  individual  properties  of  each 
stream.  The  resulting  two  combustion  disturbances  (calculated  as  volumetric  heat  release 
rate  disturbances)  were  then  simply  added  to  form  a  single  overall  disturbance.  The  overall 
disturbance  was  then  distributed  evenly  over  the  total  cross-sectional  area  of  the  augmentor, 
which  was  taken  to  consist  of  a  single  overall  stream  with  mean  mixed  properties.  This  ap¬ 
proach  accounts  for  the  different  combustion  characteristics  of  the  fan  and  core  streams, 
while  avoiding  the  complexities  associated  with  a  rigorous  treatment  of  the  radial  as  well  as 
the  axial  distribution  of  combustion  throughout  the  augmentor. 

Experience  with  modeling  the  combustion  process  as  a  plane  heat  addition  with  all 
combustion  taking  place  in  zero  length,  had  shown  that  the  resulting  predictions  of  rumble 
were  sensitive  to  the  axial  location  chosen  for  the  plane.  Since  combustion  actually  takes 
place  over  a  distance  of  30  to  60  inches,  it  was  decided  that  the  axially  distributed  nature  of 
the  burning  should  be  accounted  for.  This  was  accomplished  by  dividing  the  combustion 
zone  into  a  number  of  axial  sections,  each  of  length  t,  as  explained  in  part  1 ,  “Development 
of  Acoustic  Equations”. 

Combustion  equations  used  in  the  rumble  model  are  based  upon  an  extension  of 
empirical  steady-state  processes  to  the  case  of  time  variant  flow.  A  schematic  of  the 
steady-state  processes  is  shown  in  Figure  30.  Consider  first  that  the  augmentor  contains  only 
the  fan  stream.  An  identical  set  of  equations  will  exist  for  the  parallel  core  stream.  Following 
a  particle  of  air  as  it  moves  through  the  augmentor,  the  following  steps  will  occur: 

•  Particle  of  air  (licks  up  fuel  as  it  crosses  the  spraybar. 

•  Particle  drifts  at  through-flow  velocity  to  the  flameholder,  station  (  4  ). 

•  Particle  is  ignited  by  the  flameholder  wake  as  it  drifts  from  the  flame¬ 
holder,  to  the  beginning  of  the  combustion  zone,  station  (5)  (defined  as 
the  location  where  the  bulk  fluid  temperature  begins  to  rise  sharply ). 

•  Particle  drifts  and  burns  from  station  (5)  to  the  end  of  the  combustion 
zone,  station  (10)  (defined  as  the  location  where  bulk  fluid  temperature 
ceases  its  sharp  rise). 
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It  was  determined  (see  equation  (10))  that  for  a  linear  temperature  gradient,  the 
steady-state  volumetric  heat  release  in  a  section  of  length  1  in  the  augmentor  could  be  taken 
as  independent  of  axial  position.  This  implies  that  at  steady-state  a  particle  of  fuel-air 
mixture,  drifting  and  burning  through  the  combustion  zone,  has  a  volumetric  heat  release  rate 
that  is  independent  of  axial  position.  The  rate  can  be  computed  directly  from  the  flowrate, 
ideal  temperature  rise,  efficiency  and  combustion  zone  volume  of  the  augmentor. 

q-^WT,,  (30) 


For  small  perturbations,  it  was  assumed  that  transiently  the  volumetric  heat  release  rate  of  a 
particle  could  still  be  taken  as  independent  of  axial  position,  and  that  equation  (30)  could  be 
used  to  compute  the  rate  when  W,  T,  and  T)  are  referenced  to  instantaneous  approach  con¬ 
ditions.  The  resulting  equation  will  model  combustion  as  though  it  behaves  in  a  quasi-steady 
manner.  The  volumetric  heat  release  rate  at  any  location  in  the  combustion  zone  will  reach 
the  steady-state  value  corresponding  to  instantaneous  i  mditions  at  the  flameholder  and  at 
the  spraybar  after  a  delay.  The  delay  is  the  time  required  to  purge  the  old  combustion  gases 
and  refill  with  new  combustion  gases  traveling  at  through  flow  velocity. 

For  the  fan  stream,  instantaneous  approach  conditions  are  taken  to  be  the  instantaneous 
conditions  at  station  (3).  Because  of  the  large  pressure  drop  in  the  fuel  spraybar  injector, 
changes  in  fuel  flow  in  response  to  augmentor  pressure  at  the  spraybar  are  small  compared  to 
changes  in  air  flow.  Consequently,  fuel  flow  can  be  considered  constant,  and  the  fuel-air  ratio 
of  the  particle  as  it  crosses  the  spraybar  is  determined  by  changes  in  air  flow  only. 


.  constant 

^S/B  =  — W, — 


(31) 


A  period  of  time  =  LSC/V^  is  required  for  the  particle  to  drift  from  the  spraybar 
to  the  flameholder.  Therefore,  the  fuel-air  ratio  of  the  particle  when  it  reaches  the  flame- 
holder  can  be  expressed  as: 


FA,.(t)  -  FA.  B  (t  -  r„<  )  <32> 

At  the  ignition  plane  (flameholder)  the  particle  has  a  “potential”  volumetric  heat  release 
rate  of: 


qf  = 


£jl 

Vc 


W,  T„  Vc 


(33) 


The  ideal  temperature  rise  is  a  function  of  the  fuel-air  ratio  of  the  particlt  ffects  of  ap¬ 
proach  temperature  and  pressure  are  negligible).  The  efficiency  is  assumed  to  be  a  function  of 
the  fuel-air  ratio  and  the  approach  pressure,  temperature  and  velocity. 


T,c  -fcn(FAc) 


(34) 


Vc  =  Vfcn(FAC)  P„  T„  V.) 


(35) 


83 


The  particle  crossing  the  flameholder  will  begin  burning  after  a  time  K4/V4,  which  is  the 
time  required  to  drift  from  Station  (4)  to  (5)  while  being  ignited.  When  it  begins  burning  at 
Station  (5),  the  heat  release  rate  of  the  particle  will  be  (X  =  0  at  Station  (5)  ): 

q(o,t)  =  qr  (t  i«/V4)  (36) 

At  some  station,  X  distance  downstream  of  Station  (5),  the  local  heat  release  rate  will 
become  that  of  the  particle  after  an  additional  time  delay, 

I 

rE  =  J  dx/V(x), 


which  is  the  time  required  to  drift  from  Station  (5)  a  distance  X  at  through-flow  velocity 
V  (X).  Then,  at  a  location  X  in  the  combustion  zone,  the  heat  release  rate  will  be: 


q(x.t)  =  q(o,t  -  rE) 


(37) 


The  linearized  versions  of  equations  (31)  through  (35),  written  in  terms  of  the  Laplace 
transform  of  the  normalized  variables  are : 

FA'*  „  =  -  W; 

FA,:  =  FA  H  h  e  s 
qc  =  VV,  +  T,'  +  ri 


» -[f  ltlFA- 

:  ■[  ^  A  ]  +  [  7$],  »•  +  H  ].T-  +  [  7^1  ,V‘ 


(38) 


solving  for  qc: 

Qc  = 
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[  [fa 

jIl]  +  [ 

’  FA  J3_l  ! 

f  e  --  s 

1  L  T, 

sfaJ  c  L 

.  v  »faJ  c  J 

W, 


A  corresponding  equation  for  the  core  stream  can  be  directly  written  by  changing  sub¬ 
script  “C”  to  subscript  “H”,  and  changing  the  reference  approach  station  from  (3)  to  (3H). 


Qh 


[  >4  1 

[f*].'-  +[7wLv- 


w,„' 


(40) 
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The  total  volumetric  heat  release  rate  (subscript  “T”)  is  formed  by  adding  the  heat  release 
rates  of  the  fan  and  core  streams: 

q.vT  =  Q.  =  Qc  +  Qh  =  q.vv  +  qHvH  (41 1 

or,  in  normalized  form : 


r  q<  i 

'  Qh  1 

L  o:  Jqt  + 

-  Qt  J 

Equation  (42)  computes  the  instantaneous  volumetric  heat  release  rate  of  a  particle  of 
combined  fan  stream  and  core  stream  fuel-air  ratio  mixture  when  the  particle  reaches  the 
flameholder.  The  term  "potential”  is  applied  because  the  particle  has  not  yet  been  ignited. 
The  particle  is  ignited  by  the  flameholder  wake  as  it  drifts  a  distance  at  velocity  V^.  The 
particle  begins  releasing  the  "potential”  heat  at  station  (5),  as  defined  by  equation  (36).  To 
account  for  adding  the  core  stream  to  the  augmentor  flow  (originally  only  the  fan  stream  wa.- 
considered),  equation  (36)  was  rewritten  to  include  the  heat  release  of  both  the  core  and  fan 
streams  and  emerges  as: 

q(o.t)  =  q,  (t  — >  tf4)  (431 

Linearized : 

q  (°>t)  =  q,'  (t— #4  /  V4)  (44) 


Equation  (44)  simply  adds  a  delay  into  the  system  which  allows  tailoring  the  axial  location 
of  the  beginning  of  the  combustion  zone.  For  convenience  in  programming  the  equations, 
this  delay  is  added  to  the  drift  delay  in  the  combustion  zone  (rE)  to  form  an  overall  particle 
drift  delay  from  the  flameholder. 

=  ?«/  V4  +  rE  f 45 1 

The  particle  then  releases  heat  throughout  the  combustion  zone  as  defined  by  equation 
(37),  the  linearized  version  of  which  is: 


q'(x,t)  ~  q(o,  t  —  rE) 


1 46 1 


This  equation  represents  the  augmentor  heat  release  based  on  steady-state  conditions.  To 
accurately  model  the  combustion  process  in  the  pugmentor  the  dynamics  of  the  flameholder 
wake  must  be  included.  These  dynamics  were  incorporated  into  the  heat  release  term  as  a  first 
order  lag: 


Qiiiii/dynaniii' 


q»ut/ateadv-stalc 

1  +  tS 


The  dynamics  involved  in  calculating  the  flameholder  wake  time  constant  T  are  discussed  in 
Section  II-3. 
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Equation  (45)  was  presented  in  part  1,  “Development  of  Acoustic  Equations”,  and 
used  to  evaluate  integrals  in  equation  (7).  The  combustion  equations  require  that  the  fol¬ 
lowing  information  about  the  steady-state  operating  point  . 


■  Q<  ■ 

_  Qh  " 

’  FA  «*T.  1 

FA  ;*if 

L  Q,  J 

L  Q,  J 

L  T,  «FA  -j  c.H  • 

L 

i }  ;*FA  - 

’  P  r 

' 

- 

and 

V  .0, 

-  T)  ,lf 

C.H 

- 

n  ;>T  J  c.H 

■  n  .v 

The  heat  release  rate  ratios  Qp/Q-p  and  Q^p/Q-p  are  computed  in  the  program  from  conditions 
known  about  each  augmentor  stream : 


Qc_  =  (BFR  T„.  nc) 

Q,  (BPR  Tie  7) (•)  +  (Tl(,  TJh  ) 

Qh  _  _ <TiH  r;H) _ 

Q,  (BPR  T„  i),.)  -i  (TlH  rjH) 


(47 » 


The  partial  derivative  terms 


■  FA 

-tT 


;»■ Tj 

<7fa 


are  computed  in  the  program  from  a 

t\  H 


subroutine  curvefit  of  the  ideal  temperature  rise  curve.  A  graphical  definition  of  the  term  is 
supplied  in  Figure  31.  The  partial  derivative  terms  involving  efficiency  are  computed  in  the 
flameholder  combustion  model  and  supplied  directly  to  the  rumble  model.  Alternately,  they 
may  be  computed  from  empirical  data  and  be  input  by  the  user.  The  graphical  definition  of 
terms  is  similar  to  that  of  Figure  31. 


0.040  0.0413  0.056  0.064  0.072  0.080  0.086  0.096 

Fuel-Air  Ratio  ro  km 


Figure  HI.  Ideal  Temperature  Rise  far  Constant  Pressure  Combustion  of 
Hydrocarbon  Fuels 
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This  completes  the  combustion  equation  development  for  the  V-gutter  flameholder 
model.  All  of  the  above  equations  apply  to  the  Vorbix  and  Swirl  augmentors  except  as 
noted  below. 

For  the  Vorbix  and  Swirl  augmentors,  independent  heat  release  rates  for  the  fan  and 
core  streams  cannot  be  identified  because  of  the  flow  mixing.  In  addition,  the  effects  of 
pilot  fuel-air  ratio  on  augmentor  combustion  efficiency  must  be  accounted  for.  Equation 
(3)  is  again  applied,  but  on  an  overall  basis  only. 

Cp 

q<  =  ~  W4T,rj  (  Is  I 

The  overall  fuel-air  ratio  is  computed  from  total  mixed  air  flow  at  station  1 1  >. 

c-a  _  constant  ,  jq . 

A  '  W4 

The  overall  ideal  temperature  rise  is  a  function  of  overall  fuel-air  ratio.  The  efficiency  is  as¬ 
sumed  to  he  a  function  of  overall  fuel-air  ratio,  pilot  fuel-air  ratio  and  pressure  at  station  (  1 1 

T,  =  fch(FA)  ( 50 1 

ij  =  fch(FA,FAP,P4t 

Then  for  the  Vorbix  and  Swirl  augmentors,  the  instantaneous  “potential"  volumetric  heat 
release  rate  of  a  particle  of  mixture  when  the  particle  reaches  station  (-1 1  is: 


[  FA  aT,  1  f  FA  _a>j_  1 

q'  "  L  1  L  T,  aFA  J  l  n  oFA  J  J  V>‘ 

.  M*  ■■■■— - - 1  FAP  ,  [  liiL  1  p 
L  „  aFAP  J  L  »  ap 


Equation  (52)  applies  to  both  the  Vorbix  and  Swirl  augmentors,  and  is  equivalent  to 
equation  (42)  for  the  V-gutter  augmentor.  The  Vorbix  and  Swirl  augmentors  differ  in  pilot 
location.  The  Swirl  has  the  pilot  at  fan  duct  exit,  so  that  air  flow  through  the  Swirl  pilot  is 
proportional  to  fan  duct  exit  flow.  VV^.  The  Vorbix  has  the  pilot  near  midspan,  radially,  and 
slightly  aft  of  stations  (3)  and  (3Hl,  so  that  air  flow  through  the  Vorbix  pilot  is  proportional 
to  total  flow,  U'  |  Then,  since  fuel  flow  into  both  pilots  is  constant 

Swirl:  FAP  W, 

Vorbix:  FAP  \V4  (  ''  ' 


For  convenience  in  programming.  \V  j  can  he  replaced  by: 

vv.  w,  •  \VJH 

„.  r  bpr  ~i ...  r  i  i  ... 

v‘  I  1  •  RPR  J  4  1  •  RPR  I  XN"' 
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MICROCOPY  RESOLUTION  TESl'  CHART 

NATIONAL  BURlAll  Of  stAN[iAB[>s  1'*.-  4 


k 


Substituting  (53)  and  (54)  into  (52): 


V 

BPR  \  ( 

FAP_^_\\ 

A 

1  +  BPR/  V 

t)  aFAP  // 

-“«■{(-[“&]  -[?*])(  tUb) -(  - 

*{{'-[%&]  -[“*] )(™)N 


FA  aT, 

FA  a»i 

FAP 

aij 

V 

BPR  \\ 

L  T,  aFA  J 

ij  aFA  - 

L  *> 

aFAPJ 

A 

1  +  BPR  1) 

FA  aT, 

FA  aij 

FAP 

V 

1  )) 

L  T,  aFA  - 

if  aFA  - 

n 

aFAP- 

A 

1  +  BPR  )} 

.  P._52_lp 

L  i>  «p  J  ‘ 


Equations  (55)  and  (56)  replace  equation  (42).  All  other  combustion  equations  are 
identical  to  those  developed  for  the  V-gutter  flameholder  augmentor.  The  partial  derivatives 
in  equations  (55)  and  (56)  must  be  computed  from  empirical  data  and  be  input  by  the  user. 

This  completes  development  of  the  combustion  equations.  For  the  solution  technique, 
based  upon  applying  the  Nyquist  criterion  to  the  open  loop  transfer  function  (OLTF).  the 
OLTF  is  formed  by  renaming  q  j  to  q^  in  equation  (44)  and  by  renaming  q^  to  q()(j-p  in 
equations  (42),  (55)  and  (56). 


APPENDIX  B 

DEVELOPMENT  OF  FLAMEHOLDER  COMBUSTION  MODEL  EQUATIONS 


1.  DEVELOPMENT  OF  THE  FAN  DUCT  COMBUSTION  EQUATIONS 

The  equations  which  are  used  in  the  fan  duct  combustion  analysis  are  highlighted  in  this 
section.  The  reader  is  referred  to  the  AFWAL-TR-81-2113  (Contract  F33615-79-C-2059)  for 
full  details  of  the  analytical  development. 

The  program  utilizes  the  input  to  set-up  and  analyze  each  stream  tube  as  a  separate 
entity.  The  results  are  stored  for  final  summation  at  the  completion  of  the  fan  duct  analysis. 


The  flow  field  is  first  developed  from  the  input: 


p. 

(57) 

P*  RT. 

V.  =  M  VtRT. 

(58) 

u 

2 

ii 

£ 

(59) 

m.  =  p.V.W 

(60) 

The  streamtube  width  has  been  set  from  the  flameholder  width  and  the  blockage  ratio. 
Note  that  the  streamtube  is  assumed  to  be  1-in.  deep.  The  total  flowrates  are  thus  per  unit 
depth.  If  true  total  flowrates  are  desired,  the  number  of  streamtubes  of  each  type  must  be 
set  to  reflect  the  total  true  depth  of  that  type.  For  example,  if  5  streamtubes,  of  4  inches 
depth  each,  are  input  as  one  type,  then  set  the  input  number  of  this  type  equal  to  20. 

To  account  for  the  remova1  if  air  from  the  streamtube  for  liner  cooling,  the  input  fuel- 
air  ratio  is  adjusted  by: 


'  effective 


WCOOL 


(H 


+  BPR 
BFR 


This  increases  the  fuel-air  ratio  to  reflect  the  air  removal  when : 
WCOOL  —  mCooiinff /menKine 

BPR  =  rh,,uct/m  core 
menjine  —  nw,  "b  mcor. 

Then 

mr  =  m,  (FA)  effective 


(62) 

(63) 

(64) 

(65) 
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This  is  required  since  fuel-air  ratios  are  usually  based  on  the  total  fan  duct  air  flowrates.  If 
true  values  are  known  or  if  no  cooling  air  is  used,  set  WCOOL  =  0.0.  Refer  to  Figure  32  for 
details. 


BPR  =  WDuct/WCore  ;  WCOOL  =  WCooling/WTotal 


fO  146464 


Figure  32.  Location  of  a  Core  Streamtube  in  a  Turbofan  Engine  Augmentor 

The  injection  subroutine  divides  the  fuel  into  5  droplet  size  groups  which  represent  the 
droplet  size  vs.  volume  distribution.  The  curve  used  is  for  a  variable  area  pintle  injection. 
The  sizes  used  are: 

Group  %  Covered  Mean  Value 

1  0-20  dl0 

2  20-40  d,0 

3  40-60  d,„ 

4  60-80  d„ 

5  80-100  dK 


The  curve  is  a  function  of  the  injection  pressure  drop  where: 

APmj  =  PFSR  —  Ps  (66) 

Any  flash  vaporization  is  evaluated  from  the  fuel  enthalpy  chart  assuming  adiabatic  injection, 
i.e.,  AH  =  0 


H,  =  fen  (PFSR,  TFSR) 


(67) 


H,  =  fen  (%  vaporized,  Ps) 


(68) 
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The  droplet  vaporization  ana  acceleration  are  evaluated  by  a  small  time  step  integration 
between  the  sprayring  and  flameholder.  The  equations  are: 


dV. 


IT  ‘  tH(V- - 

for  acceleration,  and : 


(69) 


m,i„orlMd  KA,  P,  (n 


(70) 


Nu  D.  MW 
K  Rd,T. 

Nu  =  2  +  0.6  R,*'*  P,1'* 


(71) 

(72) 


for  vaporization. 

The  evaluation  of  the  liquid  temperature  follows: 


h,  =  kNu/df 

(73) 

q  =  h,  A.  (T.  -  T,)  0 

(74) 

a  _  * 

(75) 

p  e*  -  1 

z  =  Cp,  mv/irk  d ,  Nu 

(76) 

Aq  =  q  -  mv  A 

(77) 

dT(  Aq 

dt  m,  Cp< 

(78) 

4  (  d«V 

mi  =  p,~3  *  \  T  / 

(79) 

^  p.  d  f  (V.  -  V ,) 

(80) 

P« 


This  procedure  is  done  for  each  size  group  until  the  flameholder  is  reached  and  the  net 
fraction  vaporized  is  evaluated. 


0, 


•t  ?/h 


(81) 


The  impingement  of  liquid  fuel  into  the  flameholder  is  evaluated  by  use  of  a  term  02  where: 


0,  = 


m,c 


mf  | 


(82) 
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This  evaluates  the  percentage  of  the  liquid  fuel  exposed  to  the  flameholder  which  actually 
collects  into  its  surface.  The  evaluation  procedure  is  done  for  each  size  droplet  group  by  a 
correlation  of  ^  vs.  flameholder  size,  apex  angle,  flow  velocity  and  droplet  diameter.  The 
correlation  is  based  on  evaluations  performed  by  droplet  trajectory  analysis  using  the  poten¬ 
tial  flow  field  aerodynamics.  The  total  impingement  flowrate  is  thus: 

•  (83) 

A  =  JL  m(i)  d.(i) 

m  i.i 


or: 


m,c  =  d>  (1  -  di)  T  mf 


(84) 


The  liquid  film  vaporization  rate  is  evaluated  from  the  equations  for  the  surface  film  vapor¬ 
ization  caused  by  heat  transfer  from  the  flameholder  wake.  The  surface  is  broken  into  ten 
elements  and  the  vaporization  and  liquid  temperature  rise  in  each  is  calculated  from: 


m. 


C,  = 


C.  A,  P,  In  (  p^p :) 
N„  Dv  MW 


R  Ax  T. 

N„  =  0.33  R,0  *  P,°  “ 
Pv  =  fen  (T*) 


.  „  _  /  N.D.MW  \  n  .  .  /  P.  \ 

q  m,c  Cp  AT,  +  A  (  R  T>  )  P,  A,  1„  (  p>  _  py/ 


q  -  hr  A,  (Tw  TK/H) 


h.  =  if 


=  0.99  R,*  »  Pr#  “ 


(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

(91) 
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The  solution  procedure  for  breaks  the  flameholder  surface  into  10  equally  spaced  incre¬ 
ments.  The  length  of  each  is: 


Ax  = 


1  N/2 


10  sin  (a/2) 


(93) 


The  fuel  collected  by  the  surface  is  equally  divided  into  the  10  elements  on  each  face  of 
the  flameholder: 


mc(i)  =  d.  (1  ~  di)  r  m, 


(94) 
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Equations  29  to  36  are  used  for  element  i  =  1  on  the  surface  with  m  t  =  m  (i)  and  the  fuel 

*c 

temperature  is  assumed  to  be  the  same  as  the  droplet  liquid  temperature  at  the  flameholder. 
The  fraction  vaporized  is  calculated  and  the  liquid  temperature  use  evaluated.  The  procedure 
is  repeated  using  fuel  properties  evaluated  at: 


T,  (i)  =  TJ  (i)„  +-i  ATi(i)  (95) 

This  procedure  continues  until  convergence,  i.e.,  AT^  varies  less  than  1%  between  passes.  Into 
the  next  element,  i  =  2,  the  flowrate  is  set  equal  to  the  unvaporized  portion  of  the  i  =  1  flow 
and  the  collection  fraction  per  equation  (94). 


m(2)  =  mr(2)  +  mc(l)  -  mv(l) 


(96) 


This  flowrate  initial  temperature  is  set  equal  to  the  mass  average  of  the  exit  temperature 
from  i  =  1  and  the  droplet  liquid  collection  temperature: 


T„<2) 


mr(2)  T( c  +  [mc(l)  -  mv(l)l  T,t(l) 
mr(2)  +  mr(l)  -mv(l) 


(97) 


The  solution  procedure  is  separated  until  all  10  segments  are  finished.  The  vaporized  flow 
rate  is  the  sum  of  all  10  in  both  sides  of  the  flameholder: 


10 

m,  =  2  x  X]  mv(i) 
1  1 


The  fraction  vaporized,  /3;i,  is: 


mv  m  v 

m,  (1  -  0i)  r  d,  m, 


(98) 


(99) 


All  of  the  vaporized  fuel  is  assumed  to  enter  the  recirculation  zone. 

From  these  equations,  is  a  function  of  the  wake  temperature.  The  temperature  is  a 
function  of  the  wake  fuel-air  ratio  and  recirculation  rate.  Since  0$  strongly  influences  the 
wake  fuel-air  ratio,  the  solution  for  wake  composition  and  efficiency  becomes  a  curve  inter¬ 
section  procedure. 

First  we  define  the  recirculation  and  wake  kinetics  equations  and  then  the  solution 
procedure. 

a.  Recirculation 

The  wake  recirculation  flowrate  coefficient  is  defined  as: 

K,  =  mr/I' m,  (100) 

mr  =  p.  V.  N  K,  (101) 

For  mass  transfer  across  the  recirculation  zone  boundaries  and  a  homogeneous  wake: 


mr  = 


r 
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The  wake  volume  is  evaluated  as  a  function  of  blockage,  apex  angle,  and  flow  Mach  number 
from  literature  references  as  shown  in  the  flameholder  Final  Report,  AFAPL  TR-81-  .  From 
this: 


V„  -  (\  (L/DXB/DIN1 


(103) 


We  set : 


iYs 

r  '  N  (104) 

m  _  V,  P.V,, 

N  r  (105) 

Thus: 

p«V,Cv(L/D)(B/D)N 


and 

K,  =  CV(L/D)(B/D)(r )  1  ( 

By  curve  fits  of  L/D,  B/D  and  r  vs.  o,  N,  Va,  and  Ta,  we  find  the  recirculation  rate1  K  j . 

b.  Wake  Reaction  Kinetici 


The  wake  reaction  is  assumed  to  be  a  single  step,  second  order  reactive  controlled 
as  follows: 


dm  k  .  „  .  e  t:,T 

dt  R"  * °  Xf  T"  #* 


(107) 


For  a  well  stirred  reactor  (wake  is  assumed  to  behave  as  one): 


A  =  k(m-t-l)  .  ...  e~c,T 

V„P"  R>  Xo  Xf  T"  °» 


(108) 


For  the  assumed  single-step  reaction  process  postulated  here,  the  reaction  mass  balance  is 
(for  propane  fuel): 

$C,H,  +  50,  +  5mN,  —  3ycCO,  +  4ycH,0  + 

<*-y«)C,H,  +  5(1  -  y«)0,  +  5mN,  (109) 

Also,  a  linear  efficiency  vs.  temperature  function  is  assumed : 

T  =  T.  +  cAT,*,.,  (110) 


From  these  equations,  the  stinred  reactor  loading  capability  may  be  written  as : 

A  _  k(m  +  l)[5(l-y<)l*(0-y<|"~*e~C/'T'  *  'yr> 

V0P"  RByi(5(m+l)  +  <fi  +  y«l"(T,  +  <AT1"  0  ' 
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Based  on  comparison  of  predicted  results  with  available  stirred  reactor  data,  we  use  the  fol¬ 
lowing  values  for  this  reaction : 

n:  for  0<1,  n  =  20 
for  0>1,  n  =  2/0 

a:  a  =  n/2 

C:  C  =  E/R,  see  Figure  33 
This  yields: 

A  1.29  x  10”  (m  +  l)(5(l-y<)]»  (0-ye)*e  C',T<  ,115!] 

V„P*«  (0.08206)**  y«[5(m  +  l)  +  0  +  y«)**[T,  +  tAT] 


for  lean  mixtures. 


Figure  33.  Variation  in  Activation  Energy  With  Inlet  Temperature  and 
Equivalence  Ratio 

The  kinetics  solution  proceeds  by  successive  iteration  between  e  *  .999  and  0.70  to  find 
the  wake  efficiency  where: 


A  K,  P  m. 
V„P*  V„P.* 


(113) 


at  a  given  fuel-air  ratio  in  the  wake. 
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The  solution  procedure  for  the  wake  composition  and  reaction  efficiency  proceeds 
as  follows: 

( 1 )  The  wake  temperature  is  varied  in  steps  from  1000°  F  to  5000°  F  and 
calculated  at  each  wake. 

(2)  The  wake  fuel-air  ratio  is  varied  from  0.02  to  0.20  and  the  wake 
temperature  calculated  at  each  fuel-air  ratio. 

The  results  of  (1 )  are  used  in  the  wake  fuel-air  ratio  equation: 


FA -=  FA {  A  +  (1  -  <114> 

This  results  in  two  curves,  which  define  the  wake  fuel-air  ratio  vs.  wake  temperature  and 
wake  temperature  vs.  wake  fuel-air  ratio.  A  solution  technique  looks  for  the  intersection  of 
these  curves,  if  it  exists.  This  then  defines  the  stable  wake  composition  solution. 

The  fan  duct  gutter  wakes  may  be  supplied  with  hot  gases  from  an  external  (to  the  wake) 
source  such  as  a  pilot  region,  see  Figure  34.  If  this  occurs,  the  external  thermal  source  is 
assumed  to  effectively  increase  the  inlet  temperature  of  the  recirculated  air-fuel  flowrate,  i.e.,: 

K,  p.  V.  r  +  m.„  (115) 

T»  p»  v.  r  +  T,„  m„t 

m;  (116) 


mr  = 

T.  = 


Figure  34.  External  Heat  Addition  to  Fan  Duct  Gutters 
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The  program  then  analyzes  the  behavior  at  these  new  conditions  as  if  they  were  input. 

After  the  wake  has  been  analyzed,  the  turbulent  flame  penetration  into  the  free-stream 
is  analyzed. 

The  turbulent  flame  propagation  into  the  unreacted  free-stream  is  initiated  in  the 
shear  layers  of  the  wake.  The  model  used  relates  the  local  turbulent  flame  speed  to  the  loc  al 
aerothermodynamic  conditions  and  performs  a  finite  difference  integration  of  the  flame 
front  penetration  starting  in  the  wake  and  proceeding  to  the  exhaust  nozzle. 

For  the  purposes  of  current  analysis,  the  following  assumptions  were  made: 

•  Uniform  air  flow  profiles 

•  Uniform  fuel-air  ratio 

•  Incompressible  acceleration  of  free  air  velocity  by  the  fiameholder 
blockage  with  no  induced  profile 

•  Known  wake  size  and  reaction  efficiency 

•  Two-dimensional  ducted  flame. 

The  schematic  of  the  situation  which  is  analyzed  is  shown  in  Figure  35. 

The  approach  flow,  at  known  levels  of  pressure,  temperature,  velocity  and  fuel-air  ratio, 
is  accelerated  by  the  blockage  of  the  fiameholder  to  velocity  U,  where: 

u  <1  -  n 


=  Velocity  at  fiameholder  tip 
=  Approach  velocity 
=  Blockage  ratio. 

At  this  point,  Station  1,  an  induced  turbulence  level  is  calculated  from: 

■•-[{'—(  rMU]” 

This  equation  relates  the  turbulence  intensity,  eQ,  to  the  blockage  ratio  and  the  fiameholder 
zero  blockage  drag  coefficient,  C^. 


Where: 

U 

V. 

r 


9' 


Air  Velocity  Accelerated 
By  Blockage  -  Uniform 
Profile  Assumed 


xre  35.  Schematic  of  Flame  Spreading  Analysis 
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At  this  location,  the  turbulent  flame  velocity  calculations  are  initiated.  The  equation 
used  for  the  local  flame  speed  is  the  Karlovitz  equation : 

St  =  S<  +  (2d  SO* 

Where: 

St  ~  Turbulent  flame  speed,  ft/sec 
S  (  ~  Laminar  flame  speed,  ft/sec 
il  ~  RMS  turbulence  velocity,  ft/sec. 

The  value  of  u '  is: 

u'  =  «„U  (117) 

Additionally,  the  turbulent  flame  speed  initial  value  is  related  to  the  degree  of  initiation 
of  the  flame  speed  initial  value  is  related  to  the  degree  of  initiation  of  the  flame  front  by  the 
wake  by  the  following: 

St'  =  St  X  jjw  (118) 

This  generates  an  effective  turbulent  flame  speed  which  completely  fills  the  depth  of 
the  duct  and  propagates  at  the  same  transverse  rate  as  the  full  flame  speed  which  does  not 
fill  the  duct.  This  arises  from  the  fact  that  the  inefficiencies  of  the  wake  reaction  generate 
localized  regions  where  flame  front  ignition  does  not  occur.  This  use  of  a  reduced  value  ef¬ 
fective  flame  speed  accounts  for  this  in  a  two-dimensional  model. 

The  initial  value  for  the  augmentor  efficiency  is  the  wake  reaction  level  on  a  mass 
weighted  basis,  expressed  as  an  equation  this  is: 


ik 


m. 

Where: 

-  Initial  efficiency 

~  Wake  efficiency 

mr 

~  Wake  mass  flowrate 

m. 

~  Total  duct  flowrate. 

(119) 


The  type  of  flame  utilized  in  this  model  is  a  zero  thickness  flame  which  separates  a 
region  of  unreacted  propellants  from  a  region  of  completely  reacted  products.  From  this 
setup,  the  average  local  augmentor  efficiency  is  simply  the  ratio  of  the  transverse  flame 
penetration,  Ay,  to  the  duct  width,  w. 


To  be  consistent,  the  transverse  location  of  the  flame  front  at  the  initial  calculation 
station  is  taken  to  be: 


Ay„  =  T)Co  ■  w 


(120) 
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This  value  is  assigned  to  the  first  axial  station.  This  is  assumed  to  occur  halfway  down 
the  length  of  the  recirculation  zone.  From  visual  observations  of  wake  stabilized  flames, 
this  is  the  approximate  location  of  t/ansverse  flame  initiation. 

From  this  location  downstream  to  the  exhaust  nozzle,  the  flame  front  transverse  loca¬ 
tion  is  calculated  by  a  finite  difference  integration  of  the  local  flame  speed.  Several  axial 
profiles  are  introduced  as  the  integration  proceeds.  These  are: 

(1 )  The  turbulence  intensity  is  decayed  from  the  value  generated  at  the 
aft  flameholder  lip  at  a  rate  inversely  proportional  to  the  square  root 
of  axial  distance  over  an  effective  jet  length.  The  final  value  is  set  at 
the  initial  turbulence  level.  The  effective  jet  length  is  set  at  10  L/D 
where  the  D  is  the  open  area  distance  between  adjacent  flameholders. 

(2)  The  velocity  of  the  unreacted  fuel-air  mixture  is  retained  at  the  level 
generated  at  the  flameholder  lip.  Measured  profiles  from  several 
ducted  flame  test  rigs  support  this  assumption. 

(3)  A  term  is  introduced  which  relates  the  local  flame  speed  to  the  local 
average  duct  combustion  efficiency,  peaking  at  50%.  This  treats  the 
counteracting  influences  of  reduced  heat  loss  as  efficiency  increases 
and  reduces  the  free  oxygen  concentration.  Local  rates  which  roughly 
follow  a  sine  wave  function  have  been  reported  from  duct  data. 

An  additional  term  is  added  to  account  for  the  reduction  in  flame  speed  of  a  fuel  spray 
compared  to  a  premixed  flame.  This  term  relates  the  ratio  of  effective  flame  speed  to  premixed 
laminar  flame  speed.  It  accounts  for  the  complicated  interactions  during  flame  spreading  in  an 
evaporating  spray  in  a  simplified  manner.  The  effect  of  the  liquid  droplet  diameter  is  shown 
in  Figure  36.  The  droplet  diameter  utilized  in  the  analysis  will  be  the  mean  diameter  as  it 
exists  at  the  flameholder  trailing  edge. 


.  i 
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Figure  36.  Flame  Speed  for  Manodisperse  Tetralin  Spray 
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Analysis  of  the  terms  utilized  for  evaluation  of  the  laminar  flame  speed  term,  St,  has 
resulted  in  the  following: 


Where : 

Sc  =  laminar  flame  speed  at  1  atm  and  540°R 
<j>  =  equivalence  ratio 
T»  =  air  temperature,  °F 
Xoi  =  oxygen  mole  fraction. 


The  influence  of  pressure  is  indeterminate  at  this  time  and  has  been  incorporated  as 
v  Ps  for  subatmospheric  data  and  no  influence  for  pressures  above  1  atmosphere. 

The  finite  difference  solution  uses  1  in.  increments  in  axial  length  as  the  stepping  vari¬ 
able.  This  sets  a  time  interval: 

At  =  0. 08,‘13/ V.  (122) 

The  transverse  flame  penetration  distance  is  thus: 


Ay  =  St  At  =  y(i  +  l)  -  y(i) 


1123) 


where  St  is  evaluated  at  the  conditions  of  x  =  x(i). 

The  stepping  procedure  terminates  when  either  : 

(1 )  x  (i  +  1)  =  augmentor  length 

(2)  y  (i  +  1)  =  w. 

The  first  defines  nc  at  the  exhaust  nozzle,  the  second  defines  100%  rjc  before  the  nozzle. 
This  defines  one  fan  streamtube.  The  exit  temperature  is  thus: 

T„(i)  =  T„(i)  +  7r(i)  AT,(i)  (124) 

AT i(i )  =  fen  (T„(i),  FA(i)  effective)  • 

This  represents  the  actual  combustion  efficiency  based  on  the  true  fuel-air  ratio  in  the 
streamtube. 

For  multi-streamtube  cases,  the  exit  and  inlet  conditions  are  mass  averaged  using  the 
general  equation : 


The  average  input  fuel-air  ratio  and  average  inlet  temperature  combine  to  yield  tin- 
average  ideal  temperature  rise.  The  average  inlet  and  exit  temperatures  yield  the  average 
actual  temperature  use.  Thus: 


T)t.  = 


aT 

*-*  *■  HCtU 

AT  Idei 


( 1 2<; » 


This  is  the  chemical  efficiency.  The  thermal  exit  efficiency  assumes  that  the  augmentor  Icier 
cooling  air  flow  is  included  in  the  average  exit  temperature: 


L  m(i)  Tt,(i)  +  mr(10|  T„ 

f  =  Aii _  <  1 27  > 

•p.lt  n 

L  m(i)  +  nvuoi 

i  i 

This  reduces  the  average  exit  temperature  and  yields  the  lower  value  for  thermal  com- 
bustive  efficiency.  This  value  for  nc  reflects  the  average  exit  temperature  based  on  the 
average  input  fuel-air  ratio  and  based  on  total  fan  duct  air  flow  and  fuel  flow. 


Before  execution  of  the  core  streamtube  analyses,  the  influence  coefficients  which  arc 
required  are  evaluated.  These  are  of  the  form: 


—  =  Z(A)  (128) 

>'A  T) 

Where: 

A  =  V.,  p„  T„,  and  FA. 


They  are  calculated  from  a  1%  change  in  the  variables  and  the  linear  form: 

_A  _  i?a  -  7},  (A,  +A,)  (1201 

AA  r)  A2  -  A,  (tj,  +  if. 


Where : 


A,  =  1.01  A,.  (130) 

The  value  of  t?2  is  obtained  by  execution  of  the  analysis  at  all  the  same  input  as  t?  j,  except 
Aj  is  replaced  with  A2-  Thus,  the  analysis  is  done  once  for  base  and  four  more  times  for 
the  Z  factors. 
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2.  DEVELOPMENT  OF  THE  CORE  STREAM  COMBUSTION  EQUATIONS 


The  same  basic  analysis  procedure  as  accomplished  in  the  duct  is  used  in  the  core  with 
several  major  operational  differences: 

a.  There  is  no  cooling  air  removal  from  the  core  streamtubes.  Thus, 
the  input  fuel-air  ratios  are  used  in  the  analysis. 

b.  The  droplet  vaporization  rate  is  so  rapid  that  the  fuel  exists  only 

as  a  vapor  after  a  couple  of  inches  from  the  spraybar.  This  removes 
the  requirement  to  solve  for  the  wake  compositive  since  it  is  the 
same  as  the  input  fuel-air  ratio. 

c.  The  wake  reaction  efficiency  is  solved  directly  at  the  input  fuel-air 
ratio  and  recirculation  rates  which  are  calculated  the  same  as  the 
fan  duct. 

d.  There  is  no  droplet  size  effect  in  the  turbulent  flame  speed  model. 

The  rapid  droplet  vaporization  results  in  gaseous  phase  turbulent 
flame  penetration. 

The  solution  for  a  core  streamtube  proceeds  as  follows: 

(1 )  The  set-up  equations  are  the  same  as  the  fan  streamtubes. 

(2)  The  recirculation  coefficient,  Kj,  is  calculated  the  same  way  as 
done  in  the  fan  stream.  This  generates  the  value  of  A/VQP2  re¬ 
quired  for  the  kinetics  solution. 

(3)  The  wake  reaction  kinetics  solution  is  performed  at  the  same 
value  of  fuel-air  ratio  as  input  for  the  streamtube. 

(4 )  The  turbulent  flame  penetration  solution  is  the  same  as  for 
the  fan  stream  except  that  the  droplet  correction  term  is 
absent.  The  equation  introduces  a  value  for  the  oxygen 
concentration,  \D 

2 

This  value  is  less  than  the  fan  duct  due  to  the  removal  of  oxygen  by  the  mainbumer 
combustion  process.  This  vitiation  yields: 

Xo,  =  0.21  (131) 

"  *w«|olch 
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The  analysis  produces  a  value  of  rj  for  each  streamtube,  i,  by  the  same  equation  as  used  in 
the  fan: 


IJrU)  = 


Y(i) 

w(i) 


(132) 


where  Y(i)  is  the  penetration  distance  transverse  to  the  flow  and  w(i)  is  the  streamtube  width. 

The  exit  temperature  calculation  is  different  from  the  fan  duct  due  to  the  vitiation  of 
the  approach  air  flow  and  the  temperature  removal  in  the  turbine  between  the  main  com¬ 
bustor  and  the  augmentor  inlet. 


The  ideal  temperature  rise  for  each  streamtube  is  evaluated  by  generating  a  fictitious 
main  combustor  inlet  temperature.  The  procedure  is  as  follows: 

(1)  For  known  main  burner  FA  and  streamtube  inlet  temperature,  Ta(i),  a 
fictitious  AT  is  read  from  a  curve  as  in  Figure  37. 


Figure  37.  Fictitious  Temperature  Rise  vs  Main  Burner  Fuel-Air  Ratio 
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(133) 


(2)  A  fictitious  main  burner  inlet  temperature  is  calculated: 

T  mi(i)  =  T,(i)  —  ATric-t(i) 

(3)  An  overall  fuel-air  ratio  is  calculated: 

FA„.(i)  =  FA)m„  +  FA(i)  (134) 

(4)  With  FA  =  (FA)oa(i)  and  T  =  Tmg(i),  the  overall  effective 
temperature  rise  is  read  from  the  ideal  temperature  rise  curve. 

(5)  The  streamtube  exit  temperature  is: 

T,,(i)  =  AT,(i)  +  Tma<‘)  (135) 

(6)  The  streamtube  net  ideal  temperature  rise  is  thus: 

AT't(i)  =  T,,(i)  -  T.(i)  (136) 

This  value  is  calculated  for  each  streamtube  and  used  exactly  as  the  ideal  AT  cur\e  is 
used  in  the  fan  streams.  The  streamtube  exit  temperature  is: 

=  T„(i)  +  „r  ATi(i)  (137) 

The  inlet  temperatures  and  fuel-air  ratios  are  mass  averaged  as  is  the  exit  temperature, 
using  equation  (125). 

The  overall  core  efficiency  is  calculated  from  steps  (1)  to  (6)  using  average  inlet  con¬ 
ditions  to  yield  the  average  ideal  AT  and  equations  (137)  and  (125)  for  the  average  exit 
temperature: 


\T 

*  actual 

-  "^eill  T„ 

(138) 

- 

A’l'actunl 

(139) 

rji  “ 

At, 

The  influence  coefficients  are  shown  in  equations  (128)  to  (132)  are  evaluated  as  was 
done  in  the  fan. 
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